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ik  class  of  electrceagnetic  boundary  value  probless  in  which 
the  geometry  say  be  divided  into  two  seal-infinite  regions  can  be 
solved  exactly  by  the  Wiener-Kopf  technique.  Examples  of  such 
geoeetri.es  ,  called  basic  ’rfiener-Kopf  geoeetri es,  are  a  parallel-plate 
waveguide  bifurcation  and  an  infinite  array  of  thin  plates.  Certain 
of  these  basic  'rfiener-Kopf  probless  can  also  be  formulated  in  terns 
of  infinite  sets  of  linear  equations  which  say  be  solved  exactly  by 
the  conventional  residue-calculus  technique.  In  general,  however, 
the  probless  admitting  exact  solution  by  either  the  Wiener-Kopf  or 
conventional  residue  calculus  techniques  are  restricted  to  certain 
highly  idealized  geaaetries. 

A  variety  of  boundary  value  problems,  related  to  basic  Wiener- 
Hopf  geometries  but  of  sore  practical  interest,  may  be  solved  by  a 
modification  of  the  conventional  residue-calculus  technique,  by 
the  generalized  scattering  matrix  technique,  or  by  a  combination  of 
the  two  methods.  Solution  of  a  class  of  infinite  sets  of  equations 
by  the  modified  residue-calculus  technique  (MRCT)  is  based  on 
construction  of  a  meroncrphic  function  f(w )  satisfying  certain 
criteria  determined  by  the  form  of  the  equations  as  well  as  an 
auxiliary  requirement  imposed  by  the  edg.e  condition.  The,  solution 
is  obtained  by  integrating  certain  functions  related  to  f(v)  over 
contours  in  the- complex  w-plar.s. 


A  rapidly  convergent  numerical  method  for  constructing  f (w)  in  the 
MRCT  is  presented.  The  MRCT  solution  Is  given  in  a  fore  convenient 
for  computation,  not  involving  numerical  integration  or  solution  of 
large  order  matrix  equations  which  are  time  consuming  and  potentially 
inaccurate  from  a  computational  standpoint.  In  addition,  the  KRC7 
automatically  guarantees  satisfaction  of  the  edge  condition  and 
includes  built-in  tests  on  the  convergence  of  the  solution. 

An  additional  group  of  codified  Wiener-Hopf  geometries  in  which 
modal  expansions  are  possible  are  solved  by  a  coabinat  .on  of  the  MRCT 
and  the  generalized  scattering  matrix  procedure.  A-  ziattsring  matrix  . 
description,  including  both  propagating  and  evanescent  modes,  for  a 
modified  Wiener-Kopf  geometry  is  obtained  in  a  convergent  Neumann 
series  form  by  considering  the  phenomenon  of  multiple  reflections 
between  junctions  in  a  related  auxiliary  geometry.  Each  of  the 
auxiliary  junction  problems  is  solved  either  by  exact  methods  or  by 
the  MRCT. 

The  analytical  methods  described  in  this  thesis  are  applicable 
to  waveguide  discontinuities,  including  steps,  bifurcations,  and 
diaphragms,  in  both  rectangular  and  circular  waveguides;  to  phased 
arrays  of  dielectric-filled  waveguides  with  thick  walls;  to  a  variety 
of  diffraction  surfaces  and  gratings;  and  to  many  other  modified 
Wiener-Hopf  geometries.  Numerical  examples  of  discontinuity, 
array,  and  grating  problems  are  included  to  illustrate  the  application 
of  the  methods . 
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I.  INTRODUCTION 

It  is  well  known  that  a  class  of  boundary  value  problems  may  be 
solved  exactly  by  the  Wiener-Hopf  technique  (Noble  [19583).  One 
method  used  in  the  solution  of  such  problems  by  the  Wiener-Hopf 
technique  is  formulation  in  terms  of  a  semi- infinite  range  integral 
equation  of  the  form 

oo 

|  <f>(u)K(v-u)du  =  ,^(v)  (1.1) 

o 

where  6(u)  Is  the  unknown  function  and  the  kernel  function  K(v)  and 
the  inhomogeneous  term  ^iv)  are  known  from  the  geometry  and  the 
excitation.  The  solution  of  integral  equations  of  the  type  (1.1)  is 
discussed  in  the  classical  paper  by  Wiener  and  Hopf  [19313. 

An  t-i-native  procedure,  often  referred  to  in  the  literature 
as  Jones’  method  (see  Noble  [19583),  leads  directly  to  a  Fora>ier- 
transformed  version  of  (1.1)  in  the  general  form 

K(ct)  A+  (a)  =  ¥(a)  +  3.  (a),  (1.2) 

where  a  .is  the  Fourier  transform  variable.  Equation  (i.2)  has  two 
unknowns,  A+  and  B_,  analytic  in  the  upper  and  lower  half  planes, 
respectively,  of  the  complex  ct-plane  with  a  common  strip  of  overlap. 
K(a)  and  \p(.a)  are  known  functions. 

An  entirely  different  method  for  formulation  and  solution  has  been 
applied  to  a  number  of  problems  solvable  by  the  Wiener-Hopf  technique. 
The  alternative  approach  is  based  on  expansion  of  fields  in  terms  of 
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normal  modes  of  the  seal-infinite  regions  and  matching  the  resultant 
expressions  at  the  plane  of  the  interface.  Fourier  analysis  of  these 
field  expressions  leads  tc  an  infinite  set  of  simultaneous  linear 
equations  for  the  aode  coefficients.  The  exact  solution  of  these 
equations  by  a  residue-calculus  technique  is  discussed  in  Section  II. 
An  alternative  method  of  solution  by  direct  inversion  has  been  found 
by  Mittra  C1S63}. 

The  Wiener-Hopf  method,  and  the  residue-calculus  technique  are 
quite  versatile,  finding  application  is  father  diverse  areas  of 
mathematical  physics,  e.g. ,  acoustic  and  electromagnetic  wave 
propagation,  design  of  optimal  filter's,  etc.  Also,  the  Wiener-Hopf 
and  residue  calculus  techniques  are  the  only  known  methods  for  exact 
solution  of  a  class  of  problems  not  conforming  to  separable  coordinate 
systems.  Nonetheless,  the  scope  of  these  techniques  is  United  to 
idealized  geometries  satisfying,  certain  strict  requirements . 

In  recent  years,  there  have  been  attempts  toward  extending 
Wiener-Hopf  and  residue-calculus  techniques  to  a  wider  class  of 
problems  of  more  practical  interest.  Thess  may  be  described  as 
modified  Wiener-Hopf  geometries,  since  they  can  be  related  to  a 
basic  Wiener-Hopf  geometry.  Typical  examples  are  diffraction  by  a 
thick  half -plane,  scattering  at  the  junction  of  two  cylindrical  or 
parallel-plate  waveguides  of  dissimilar  transverse  dimension, 
radiation  from  a  dielectric  loaded  waveguide  phased  array,  etc.  In 
this  thesis,  some  methods  are  presented  for  formulating  and  solving 
a  class  of  problems  of  the  modified  Wiener-Hopf  type. 


A  rapidly-copvergent  numerical  method  based  upon  an  extension  of 
the  conventional  residue-calculus  technique  is  presented  for  the 
solution  of  a  class  of  infinite  sets  cf  simultaneous  linear  equations. 
In  this  method,  the  need  for  numerical  inversion  of  large  matrix 
equations,  such  as  those  derived  in  an  integral  equation  approach  to 
these  boundary  value  problems,  is  eliminated.  In  the  present 
approach,  the  solution  is  obtained  in  the  same  form  as  generated  by 
the  conventional  residue-calculus  method,  permitting  rapid  numerical 
calculation.  In  addition,  the  procedure  guarantees  satisfaction  of 
the  edge  condition,  which  is  a  significant  advantage  over  certain 
iterative  techniques  where  it  is  difficult  to  advance  a  direct  proof 
that  the  edge  condition  is  indeed  satisfied. 

A  class  of  related  problems  is  solved  by  the  generalized 
scattering  matrix  procedure  (Hittra  and  Pace  [1963]),  in  which  the 
method  of  multiple  reflection  between  junctions  is  extended  to 
inclr  le  some  junctions  which  are  not  exactly  solvable,  but  for  which 
very  accurate  numerical  solutions  can  be  obtained  by  the  modified 
residue  calculus  technique.  For  the  sake  of  completeness,  a  brief 
description  of  the  scattering  matrix  procedure  is  included  in  Section 
II.  Applications  to  phased  arrays  as  well  as  diffracting  surfaces 
and  gratings  are  made  in  Sections  IV  and  V,  illustrating  that  a 
combination  of  the  modified  residue-calculus  method  and  scattering 
matrix  techniques  is  useful  for  resolving  a  rather  large  variety  of 
boundary  value  problems  which  may  be  identified  as  modifications  of 
the  basic  Wiener-Hopf  geometry. 
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II.  ANALYTICAL  METHODS 


The  analytical  methods  applied  t-b  the  problems  in  this  thesis  are 
described  in  this  section.  The  method-  ar-e  applicable  to  3  class  of 
modified  Wieuer-Hopf  type  geometries  in  which  the  electromagnetic 
fields  may  be  expanded  in  terms  of  a  discrete  set  of  modes.  The 
eigenvalue  spectrum  associated  with  these  problems  is  discrete  by 
virtue  of  either  restriction  to  closed  regions  of  periodicity  i-i  open 
regions.  In  this  event,  the  fields  are  completely  described  by  the 
set  of  mode  coefficients.  The  methods  described  here  are  used  to 
obtain  the  required  mode  coefficients  by  a  rapidly  convergent 
numerical  technique. 

The  Modified  Residue-Calculus  Technique 

The  modified  residue-calculus  technique  (MRCT)  recently  introduced 
by  Mittra,  Lee,  and  VanBlaricum  [1968],  is  a  method  for  the  approximate 
solution  of  a  class  of  infinite  sets  of  simultaneous  linear  equations, 
such  as  those  generated  by  a  mode  matching  procedure  at  a  discontinuity 
in  a  waveguide.  When  the  equations  admit  an  exact  solution,  the  MRCT 
gives  that  solution,  but  the  approximate  solutions  to  other  problems, 
is  obtained  without  the  need  for  numerical  inversion  of  extremely 
large-order  matrix  equations. 

The  appropriate  sets  of  equations  have  the  general  form 


£  A  (— = — 
n=I  " 


)  =  (  1 


m 


a  +y 
n  'm 


a  +y 

p  m 


a  -y 
p  m 


•),  m  =  1,  2,  3,  ...  (2.1) 


where  {a  },  {y  k  {A  }  are  all  known  and  {A  }  are  to  be  determined.  The 
n  n  n  n 

MRCT  is  an  extension  of  the  conventional  residue-calculus  technique,  which 


- 4. 


-  ^ 


provides  an  exact  solution  to  a  set  of  equations  of  the  general  fora 
(2.1)  when  A  =  0,  o  =  1,  2,  3,  .... 

D 

The  conventional  residue-calculus  technique  was  introduced  by 
Berz  [13513  and  applied  to  the  set  of  equations  occurring  in  the 
problem  of  diffraction  of  a  plane  wave  by  a  set  of  parallel  plates. 
Subsequent  applications  of  the  technique  have  been  cade  by  Whitehead 
[19513,  Hurd  and  Gruenberg  [19543,  and  Wu  and  Galindo  [19563  "to 
cases  which  are  relevant  to  certain  aspects  of  this  thesis.  Another 
application  of  the  method  is  found  in  Collin  [19603. 

The  conventional  residue-calculus  technique  can  perhaps  best  be 
described  by  an  example  which  is  related  to  a  set  of  -  _itions  derived 
in  Section  III.  These  equations,  related  to  a  sit pie  waveguide 
bifurcation,  are 


Z  Vo  -6  }  "  a  +6  ,  m  =  1,  2,  3,  ... 
n=l  n  m  pm 


(2.2) 


E  A  (-— 1--  )  =  —4 —  ,  m  =  1,  2,  3,  .. 
,  n  a  -y  a  +y  *  ’  ’  * 

n=l  n  m  p  'm 


(2.3) 


where 


<  2-(n7i/a)2  =  -  j  /(nu/a)2-k  2 
o  o 


B  =  rk  2-(nti/b)2  =  -  j  *4nu/b)2-k  2 
no  o 


y  =  /k  2-(nt/c)2  =  -  j  /(nn/c)2-k  2 
no  o 


and  the  equations  are  to  be  satisfied  simultaneously. 


It  is  desired  to  construct,  by  means  of  contour  integrals  of  a 


known  meromorp'nic  function,  equations  of  exactly  the  same  form  as  (2.2) 
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and  ( 2 r3 )  and  to  solve  for  the  {A^}  by  Identifying  corresponding  terns. 
The  required  nerooorphic  function  g(w)  satisfies: 

Cg:l)  g(w)  has  sinple  poles  at  v  =  ®0,  «3,  ...»  and  at 

w  =  -  a  . 
p 

(g:2)  g(Bn)  =  0,  n  =  1,  2,  3 . 

(g:3)  8(Yn)  =  0,  n  =  1,  2,  3,  .... 

-3/2 

(g:4)  g (w)  is  algebraic  (specifically  g(w)  ^  Kw  )  as 

jw|  ■+  - 

Condition  (g:4)  is  a  physical,  constraint  of  the  problen  (in  this  case, 
an  edge  condition)  and  insures  a  unique  solution.  Techniques  for 
constructing  such  a  function  are  well  known,  and  giw)  takes  the  fora 

g(w)  =  K  e*p(tw)  (2.01 

P 


where 

CD 

H(w,B)  =  n  (1  -  —)  exp(jwb/nir)  (2.5) 

n=l  n 

and  similarly  for  JI(w,cx)  with  b  replaced  by  a  and  JI(w,y)  with  b 
replaced  by  c.  Also, 


L  =  [£n(a/b)  +  ~  £n(b/c)] 


(2.6) 


In  the  complex  w-plane,  let  C  be  a  circle  of  infinite  radius 
enclosing  all  of  the  poles  and  zeros  of  g(w),  and  consider  contour 
integrals  of  the  type 


JL  l  sL 

2nj  J  w- 


g(w) 


dw ,  m  =  1 ,  2,  3, 


m 


and 


8 


As  before,  construct  a  mercaorphic  function  f(w)  having  specified 
properties: 

(f:l)  f(w)  has  simple  poles  at  w  =  a, ,  a^,  a^,  ...»  and  at 

w  =  -  a  . 

P 

(f:2)  f(8n)  =  0,  n  =  1,  2,  3,  .... 

(f:3)  f(yn)  +  '.nf(-Yn)  =  0,  r.  =  1,  2,  3,  .... 

(f:4)  f(w)  is  algebraic  as  |w|  -*■  ®.  Specifically,  f(w)  ^  Kw 
1  <  v,  and  the  exact  value  of  v  is  determined  by  a 
physical  constraint.  It  will  be  shown  later  that  v  may 
be  determined  directly  from  the  asymptotic  behavior  of 
the  X  *s  as  m  becomes  large. 

TSi  . 

It  will  be  noted  that  conditions  (f:l)  and  (f:2)  are  identical 

to  (g:l)  and  (g:2)  imposed  on  g(w),  the  function  appropriate  for  the 

reduced  problem  with  A  =0.  Condition  (f:3)  is  a  generalization  of 

m 

(g:3). 

If  such  a  function  f(w)  satisfying  (f :1~4)  can  be  constructed, 
then  the  solution  to  the  equations  proceeds  as  before.  Consider  the 
contour  integrals  along  C 
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Z  R  (a  )  (— — 
,  f  n  a  -y 
n=l  n  m 


a +Y*  '  Rf("ap)  (a+Y, 


-  + 


a 


n  m 


p  m 


a  -y 
p  m 


♦  +  V!-Tn]  *  °*  °  *  x» 2’ 3' 


(2.13) 


*  W  (rV)  -¥"%>  (r?r>  + 

n=l  n  m  *  pm 


f(8) 

o 


=Otn=l,  2,  3, 


(2.14) 


where  R,_(a  )  =  residue  of  f(w)  at  w  =  a  .  Provided  f(v)  is  normalized 
r  n  n 

to  make  R,(-a^)  =  1,  (2113)  and  (2.14)  may  be  identified  with  (2.11) 
and  (2.12)  yielding 


An  =  W’  n  "  1>  2>  3»  •••  (2,15) 

Construction  of  f(w).  The  key  step  in  the  modified  residue- 
calculus  technique  is  the  construction  of  f(w).  Unlike  in  the 
conventional  technique  where  all  of  the  poles  and  zeros  of  the  function 
as  well  as  its  asymptotic  behavior  are  prescribed,  in  the  modified 
technique  only  a  subset  of  the  collection  of  zeros  is  known  explicitly, 
while  all  of  the  poles  and  the  asymptotic  behavior  are  again  known. 
Since  f(w)  must  coincide  with  g(w)  in  the  limiting  case  =  0,  f(w) 
may  be  written  in  the  following  form: 


«»>  *  *-•«*>  • 


(2.16) 


where  L,  n(w,8),  and  n(w,a)  are  given  by  (2.5)  and  (2i6),  and  where 


H(w,y')  =  II  (1  -  7^-p)  exp(+jwc/niT), 
n=l  fn 


[2.17) 


with  {yn’ }  being  the  as  yet  unknown  zeros  of  f(w).  Comparison  of 
f(w)  with  g(w)  shows  that  except  for  constant  factors  the  two  functions 
are  related  by  shifting  an  infinite  set  of  zeros  (Yn)  of  g(w)  to 
{Yn*}  of  f(w),  and  that  when  X^  =  0,  the  zeros  must  coincide. 

The  asymptotic  values  of  the  set  {■*  ’}  may  be  determined  by 
imposing  condition  (f:3).  Let  =  Yn’-Yn  be  the  difference  of  the 
zeros.  Condition  (f:3)  requires: 

JI(y  ,S)JI(y  .Y')  K'  exp(-LY_)n(-YTn»S)n(-Y  ,y') 

vi  /-  .  m  m  ,  ro  m  m 

*Xp  Ym  (y  +a  )H(y  ,o)  (-y  +a  ) II ( —  y  »a) 

ro  p  m  m  p  m 

=  0,  n  =  1,  2,  3,  ...  (2.18) 

For  m  large,  the  constants  may  be  replaced  by  their  asymptotic  values, 
i.e.,  &n  -jnir/b,  ^  -jmr/a,  Yn  ^  -jnm/c,  and  ^  D  (an  assumed 
asymptotic  value),  so  that  (2.18)  becomes  approximately 


K.  .  CC1  -  ^exp(S|)]  «  CU  -  Tgfojh*®! 

n=l  n=l 


/  •  •  \  TT  ?"  /  T  nic  \  -  /fllG  - 

(-]—  +  a  )  n  L (1  -  — )exp(— )] 
c  p  ,  na  na/J 

*  n=l 


+  X  K' 
m 


exp(+jL^)n  [(1  +  g)exp(- g) 
n=l 


C(1  t  S» 


=  0,  m  =  1,  2,  3,  ....  (2.19) 

Equation  (2,13)  can  be:  reduced  by  observing  that  the  infinite 


products  take  the  form  of  an  infinite  product  expansion  for  the  gamma 
function  (see  Noble  [1958]). 


L2  sin(-xjr/a-b7r/a)  exp(-Cx/a)T(H-b/a) 

r(l+b/a-x/a) 

for  x  -  -  |x|  (2.20b) 

where  C  =  Euler's  constant. 

Utilising  these  relations  in  (2.19),  one  has  after  rearranging  and 
canceling  like  terms: 

2  sin(-imr-jcD)  =  +  X  sin(-mna/c)/sin(-nnrb/c)  (2.21) 

m 

But  sin(-nnra/c)  =  sin(-mir-mtb/c)  -  (-l)m  sin(-imrb/c) 

and  sin(-mTr-jcD)  =  (-1)  sin(jcD),  so  that  (2.21)  becomes,  in 

the  limit  as  ra  +  », 


2  sin(jcD)  =  -  lim  X 

m 


or  D  =  ■—  sin_1(-  ^  lim  X  )  (2.22) 

m-*° 

The  explicit  evaluation  of  D  =  lim  D  permits  a  priori  determination 

n-*° 

of  how  the  zeros  of  f(w)  are  to  be  shifted  for  large  indices.  It  will 
be  shown  subsequently  that  an  expression  equivalent  to  (2.22)  may  be 
obtained  directly  from  the  edge  condition. 
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Exact  determination  of  the  shifted  zeros  (Yn’)  for*  every  n  is  not 
possible.  However,  for  purposes  of  numerical  computation  it  suffices 
to  use  the  determined  asymptotic  value  for  all  but  the  first  few,  say 
H  of  them,  and  to  find  those  first  M  values  explicitly.  Also,  except 
for  exponential  factors  and  an  irrelevant- constant  factor  X",  the 
first  M  terms  of  the  infinite  product  may  in  general  be  replaced  by 
a  polynomial  of  degree  M.  That  is, 


M 

K* 1  II  (l-w/y^  ')  =  w^  +  a^w**  ^  +  ...  +  a^  ^  w  +  a^ 

n=l 


*  P(w) 


(2.23) 


where  to  specify  the  functional  behavior  it  is  sufficient  to  determine 
the  M  polynomial  coefficients  a^,  a2,  ...,  a^. 

Accordingly,  the  function  may  be  written 


f (w)  a  K’P(w)  exp(Lw)n(w,g)lT^M^(w,Yl ) 

(w+a  )II(w,a) 

P 


(2.24) 


Cm) 

where  Jr  y(w,Y*)  denotes  the  omission  of  the  first  M  factors  from  (2.17) 
and  replacement  of  Yn*  by  Yn  +  D  for  all  n  >  M.  For  simplification 
Write 

f(w)  =  K'  P(w)f'(w),  (2.25) 


where  f’(w)  incorporates  all  of  the  known  parts  of  f(w). 

The  polynomial  coefficients  may  be  determined  by  imposing  condition 
(f : 3 )  for  n  a  l,  2,  3,  M.  Specifically, 


r  H  A  M-l 
*  [Yn  +  alYn  + 


*  aM-l  Yn  + 


*  •  • 
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t  An  K'  C(-rn)H  +  +  ...  +  a1I_1(-Tn)  +  =  ° 

for  n  =  1,  2,  3,  ...,  M,  (2.26) 

yielding  a  set  of  M  simultaneous  linear  equations  for  a^,  a^,  . .., 

a^.  Solution  of  the  set  of  equations  may  be  accomplished  by  usual 

numerical  methods.  Once  P(w)  is  known,  the  constant  K’  may  be 

chosen  to  satisfy  R,.(-a  )  =  1. 

f  P 

It  must  be  noted,  however,  that  even  though  the  MRCT  requires 
numerical  solution  of  a  set  of  linear  equations,  the  difficulty  of 
that  solution  in  no  way  approaches  that  of  attempting  a  direct 
inversion  of  a  truncated  form  of  the  original  infinite  set.  For 
most  instances,  it  suffices  to  solve  a  5  x  5  or  smaller  set  of 
equations  explicitly  in  the  MRCT,  while  much  larger  ones  wculd  be 
appropriate  for  the  truncated  set.  In  addition,  the  truncated  set 
for  (2.11)  and  (2.12)  is  quite  ill-conditioned,  while  the  set  for 
the  MRCT  is  well-behaved.  An  example  in  Section  III  illustrates  the 
rapid  convergence  of  the  solution  via  the  MRCT  for  a  waveguide  step 
discontinuity. 

An  additional  feature  of  the  MRCT  for  numerical  work  is  that 
a  built-in  convergence  test  is  available,  because  the  function  found 
numerically  must  satisfy  conditions  ( f : 2 )  and  (f:3).  A  numerical 
test  of  these  conditions  may  be  almost  trivially  incorporated  into 
a  computer  program. 

Perhaps  the  greatest  advantage  of  the  MRCT  is  that  satisfaction 
of  the  edge -condition  via  (f:4)  is  a  priori  guaranteed.  And,  more 
importantly,  when  the  slight  modifications  of  the  geometry  alter 
the  edge  behavior,  the  correction  is  made  automatically  in  the  MRCT. 

A  further  discussion 'of  this  property  is  giver,  in  Section  III. 
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Generalized  Scattering  Xatrix  Technique 

The  generalized  scattering  matrix  Technique  is  another  method  for 
solving  a  class  of  boundary  value  problems .  The  concept  of  a  generalized 
scattering  matrix,  introduced  by  Mittra  and  Pace  [19533  is  very  closely 
related  to  the  scattering  matrix  of  circuit  theory  (Seshu  and  3alabanian 
[19593)  or  of  microwave  network  theory  (Collin  [19563).  However,  it 
differs  in  that  it  is  extended  to  consider  evanescent  as  well  as 
propagating  codes,  so  that  in  waveguide  junction  problems,  the 
scattering  matrix  will  in  general  be  of  infinite  order.  (In  subsequent 
discussions,  the  term  "scattering  matrix"  will  mean  "generalized 
scattering  matrix . "  ) 

The  scattering  matrix  is  defined  for  a  junction  between  two 

regions  in  which  the  fields  cay  be  expanded  in  modes,  such  as  the 

waveguide  bifurcation  of  Figure  2.1.  If  the  nth  7E  mode  is  incident 

upon  the  plane  z=0  from  region  (A),  fields  are  reflected  back  into  (A) 

and  transmitted  into  (B)  and  (C).  The  component  of  the  field  is 

Til 

expanded  in  a  Fourier  sice  series,  and  the  coefficient  cf  the  m 

node  referred  to  the  plane  z=0  is  defined  as  the  amplitude  of  that 

mode.  If  the  amplitude  of  the  n**1  mode  incident  in  (A)  is  one,  then 

Th 

the  amplitude  of  the  n  scattered  code  in  (A)  is  S'  (m,n)  and  the 

amplitude  of  the  m**1  transmitted  mode  in  (B)  is  SBA(a,n).  Terms 

AA  BA  AA  BA 

S  (m,n)  and  S  (m,n)  are  the  general  matrix  elements  of  S  *"  and  S  , 

respectively.  The  other  scattering  matrices  are  similarly  defined. 

In  conventional  scattering  matrix  formulations,  the  modes  are 

normalized  so  that  a  propagating  mode  carries  unit  power.  Since  the 

generalized  scattering  matrix  includes  evanescent  modes,  such  a 


normalization  is  inappropriate.  One  consequence  is  that  the 
scattering  matrices  are  not  symmetric . 


The  generalized  scattering  matrix  technique  is  a  means  of 
solution  for  problems  which  are  related  to  an  exactly  solvable 
auxiliary  problem.  For  example,  the  inhomogeneously-filied  waveguide 
bifurcation  of  Figure  2.2  is  related  to  the  simple  bifurcation  of 
Figure  2.1,  which  may  be  solved  exactly  by  either  Wiener-Hopf  or 
residue-calculus  techniques  (Section  III).  If  the  dielectric  region 
(C*)  is  recessed  as  in  Figure  2.3,  the  resultant  problem  involves 
two  junctions,  namely  junction  1,  the  original  bifurcation,  and  junction 
2.,  the  interface  between  regions  (C)  and  (C 1 ) .  Each  junction  problem 


may  be  solved  individually  and  represented  by  scattering  matrices, 
*  „  QAC  ab  -CA  „CC  CC' 

Suppose  now  that  a  TE  wave  is  incident  upon  junction  1  from 


region  (A).  This  wave  is  completely  determined  by  the  coefficients 


of  the  modes,  and  hence  can  be  represented  by  a  mode  column  vector 

"til 

<{>  such  that  is  the  coefficient  of  the  icn  mode  in  the  incident 
field.  At  junction  1,  fields  are  reflected  back  into  region  lA)  and 


transmitted  into  regions  (B)  and  (C).  The  mode  vector  for  the  field 

AA 

reflected  back  into  (A)  is  S^  <J>,  and  for  the  field  transmitted  into 

CA 

(C),  the  vector  is  S^  $.  Part  of  the  wave  transmitted  into  (C)  is 

reflected  back  into  (C)  and  part  is  transmitted  into  (C)  at  junction 

C  'C  CA 

2.  The  mode  vector  in  region  (C * )  is  S^  S^‘  <{>.  Figure  2.4  depicts 
symbolically  the  multiple  reflection  phenomena  between  junctions  1 


and  2.  The  fields  in  each  region  are  the  sums  of  the  fields 


transmitted  or  reflected  successively.  It  should  be  noted,  though, 


JUNCTION  2 


for  singly- inhomogeneous  bifurcation 


Multiple  reflection  diagram  for  waveguide  bifurcation 
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that  since  the  effect  of  propagation  over  the  distance  5  has  been 
neglected,  the  j  -suits  are  valid  only  as  6  -*■  0,  in  which  case  the 
geometry  reduces  to  that  of  Figure  2.2,  the  desired  geometry.  The 
fields  in  the  dielectric  region  may  be  represented  by  the  mode 
vector  tj>,  which  is  found  by  summing  the  fields  transmitted  into  (C*). 


-  sc2'csc^  + 


C'CCCCCCA 

S2  s1  S2  SjL  $ 


+  SC‘C(SCCSCC)2SCA^  +  #< 


=  £  SC'C(sfsf)nSCk$ 

n=0 


(2.27) 


Equation  (2.27)  takes  the  form  of  a  Neumann  series,  which  can  be 
summed  in  the  standard  way.  Proof  of  the  convergence  is  given  in 
Pace  [1964].  Also,  the  mode  vector  ip  must  be  given  in  terms  of  the 
composite  junction  scattering  matrix  (no  subscript),  so  that 


.  _  cC'Aa 
ip  =  S  <p 


(2.28) 


where  I  is  the  identity  matrix.  Since  the  result  (2.28)  is  valid 
for  arbitrary  incident  mode  vector  <p,  the  implication  is  that 


,C‘A  _  cC'CrT  _CC_CC.-1QCA 
3  S2  (I"S1  S2  '  S1 


(2.29) 


The  other  scattering  matrices  for  the  problem  can  be  found  in  an 


analogous  manner. 
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Equation  (2.29)  is  formally  exact,  since  all  of  the  auxiliary 
scattering  matrices  may  be  found  exactly.  For  practical  purposes, 
however,  results  are  found  approximately  by  truncating  the  scattering 
matrices  and  performing  the  inversion  numerically.  Many  interesting 
problems  may  be  solved  approximately  by  the  generalized  scattering 
matrix  procedure  even  when  the  auxiliary  problem  cannot  be  solved 
exactly.  Often  a  very  good  approximate  scattering  matrix  represent¬ 
ation  of  the  auxiliary  problem  can  be  found  (for  example,  by  the 
modified  residue-calculus  technique)  so  that  the  composite  scattering 
matrices  may  be  found  as  before.  An  example  of  this  kind  is  given 


in  Section  III. 
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III.  INHOMOGENEOUSLY-FILLED  WAVEGUIDE  BIFURCATION 

Among  the  problems  occurring  in  micr-owave  theory  are  waveguide 
discontinuity  problems.  Host  of  the  problems  of  interest  do  not 
possess  exact  solutions,  so  that  many  approximation  techniques  for 
the  solution  of  such  problems  have  been  developed.  Many  such 
techniques  and  their  application  to  finding  equivalent  transmission 
line  and  lumped  circuit  element  representations  for  waveguide 
discontinuities  are  given  by  Marcuvitz  [1951].  Since  usually  the 
equivalences  are  treated  with  respect  to  one  or  two  propagating 
modes,  description  of  discontinuities  is  often  given  in  terms  of 
modal  reflection  and  transmission  coefficients  (i.e.,  scattering 
matrix  elements). 

The  simple  bifurcation  problem  (Figure  2.1)  admits  an  exact 
solution,  but  the  generalization  to  the  inhomogeneously-filled 
bifurcation  (Figure  3.1),  which  may  be  applied  to  waveguide  junctions, 
waveguide  diaphragms,  and  even  to  phased  arrays,  cannot  be  solved 
analytically  in  closed  form.  However,  the  modified  residue  calculus 
technique  is  expressly  suited  for  solving  such  problems  approximately 
and.  will  be  applied  in  this  section  to  the  solution  of  the  inhomogen¬ 
eously-filled  waveguide  bifurcation.  Numerical  results  for  a  set 
of  special  cases  are  given,  and  application  to  related  problems 
are  discussed.  Application  of  the  results  to  the  phased  array  is 
made  in  Section  IV. 

Analysis 

The  geometry  of  the  inhomogeneously-filled  bifurcation  is  given 
in  Figure  3.1.  The  septum  at  x  =  c,  z  >_  0,  is  assumed  to  be  infinitely 


thin  and  perfectly  conducting.  Regions  (B)  and  (C)  are  filled  with 


media  of  relative  dielectric  constants  and  ec>  respectively. 
Assume  that  TE^q  modes  with  electric  vector  parallel  to  the  edge 
of  the  septum  are  incident  upon  the  plane  z=0  from  regions  (A)  and 
(B)  with  amplitudes  A  and  B,  respectively.  By  linearity,  the 
excitations  may  be  treated  either  together  or  separately.  No  loss 
of  generality  results  from  not  considering  explicitly  excitation 
from  (C),  since  the  role  of  (B)  and  (C)  may  be  interchanged  and 
the  principle  of  superposition  applied. 

The  total  electromagnetic  fields  may  be  derived  from  the  only 
non-zero  component  of  electric  field,  E  . 


(3.1a) 


H  =  -  n-i-  —  (E  ) 
z  jwp  3x  y 


(3.1b) 


(exp(jut)  time  variation  is  assumed  throughout).  The  analysis  and 
exact  solution  for  the  reduced  problem  -  1  was  given  by 

Hurd  and  Gruenberg  [1954].  Formulation  for  the  generalized  problem 
proceeds  in  a  similar  manner. 

The  formulation  of  the  problem  is  facilitated  by  an  auxiliary 
geometry,  that  of  Figure  3.2.  Dielectrics  in  regions  (B)  and  (C) 
are  recessed  small  distances  u  and  v,  forming  regions  (B1)  and  (3), 
and  (C)  and  (C),  respectively.  In  the  five  regions,  the  transverse 
electric  field  may  be  written 


Figure  3,2  Auxiliary  geometry  for  inhomogeneously-filied  bifurcation 
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(A)  S 

y 

O')  2y 

(3>  Ey 

(C*)  E 


Z 

n=l 


Z 

n=l 


A^  sintcxx/aJexpCJa^z)  +  A  sin(pTx/a)exp(~3o^±)  (3.2a) 


sin(u-(x-c)/b)  [3  exo(-jS  z)  +  D  exp(j3  z)3  (3.2b) 

71  *  H  ji  a 


F  sin(nx(x-c)/b)exp(-jlf  z)  +  3  sin(px(x-c)/b)exp(jg-,z) 
n  K  *  — 


(3.2c) 

sin(nsx/c)  CC  e3o(-jy  z)  +  S  exp(jy  z)3 
n  n  n  *  n 

(3. 2d) 

G  sin(nsx/c)exo(-jY  z) 
n  *  n 

(3.2e) 

where 


o  =  *4:  2*-(ns/a)2 
n  o 

3  =  *4: 2-(iw/b)^ 

n  o 

Yn  =  42-(^/c)2 

Bn  =  ^o2-(iw/b)2 

T  =  >4  k  2-(nz/c)2 

n  co 


2 


u 


U 


e 

o  o 


Except  for  a  constant  factor,  the  transverse  components  of  the 
magnetic  field  in  the  corresponding  regions  are  found  by  partial 
differentiation  with  respect  to  z. 

Hatching  tangential  field  components  at  z  =  v,  0<_x<^c,  and 
equating  Fourier  coefficients  yields 


(3.3a) 


For  u  =  v  =  0,  the  auxiliary  and  original  geometry  coincide,  so  that 
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the  relations  of  (3.5a-b)  nay  be  used  in  formulation  for  regions  (A), 
(3),  and  (C). 

Matching  tangential  field  coaponents  at  the  plane  2=0  gives 
the  following  expressions: 


A  sin  (psx/a)  +  Z  A  sin(n~x/a)  = 
n=l  n 


Z  C  (1+A  )sin(nsx/c),  0  <  x  <  c 
n=l  n  n  “  ~ 


Z  B  (1+5  )sxn(n^(x-c)/b) 
,  n  n 
n=l 


(3.6a) 


+  2 - — B  sin{pt(x-c)/b),  c  x  a  (3.6b) 

(B  +B  ) 

P  P 


-c  A  sin(p7rx/a)  +  Z  a  A  san(nitx/a) 
P  .  n  n 

r  n=l 


Z  y  C  (-1+A  )sin(.mtx/c)  0  <  x  <  c  (3.7a) 

.  n  n  n  —  — 

n=l 


Z  B  B  (-1+5  )sin(nir(x-c)/b) 
n=l  n  11  n 


+  2 — B  sin(p7t(x-c)/b)  c  f.  x  ,1  a  (3.7b) 

(0  +3) 

p  p 


Relations  among  the  mode  coefficients  may  be  obtained  by  multiplying 
the  equations  above  by  sin(mtx/c)  and  sin(mm(x-c)/b)  and  integrating 
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in  the  respective  regions  0  x  c  and  c  <_  x  <_  a.  Employing  the 
integral  identities 


sin(mux/c)sin(mrx/a)dx  =  (-l)m+^  — 

c  2  2 


a  -y 
n  m 


and 


sin(nm(x-c)/b)sin(mrx/a)dx  =  - 


.mu  'sin(mrc/a 

b  2  q2 
o  -B 
n  m 


the  relations  for  the  mode  coefficients  become 


^sin(pttc/a)  E  A  s^n^n7TC/a^ 


2  2 
a  -y 
p  m 


,  n  2  2 

n=l  a  -y 
n  m 


=  (-l)ra+1  c  (1+X  ) 
2mn  m  m 


_  Asin(pirc/a)  „  _  A  sin(n7rc/a) 

V  2  va  +  E,  “n  An  2  2 

a  -y  n=l  a  -y 

p  m  n  'm 


=  (-l)m+1  ~~  Y  C  (-1+X  ) 
2nnr  mm  m 


Asin(pKc/a) 

A  2  2 
a  -B 
p  m 


+  Z  A 


sin(niTc/a) 


.  n  2  „2 
n=l  a  -B 
n  m 


B  (1+5  ) 
2mir  m  m 


,2  2gr 


0  - B6P 

2p"  (6  +S  )  " 

P  P 


(3.8) 


(3.9) 


(3.10) 


I 
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t  » 

tv 


Aslnfec/a)  £  a  A 

P  <£-S*  n=l  »  "  aUl 

pro  n  ro 


~  &  B(-l+5) 
zmir  ro  in  ro 


,2  23  3 
b - E_P-B6P 

2p71  (3  +3  )  “ 

P  P 

(3.11) 


For  simplification  in  notation,  let 


B  =  2 - 2—  B 

(3  +3  ) 

P  P 


A’  =  A  sin(nnc/a) 
n  n 


(3.12a) 


(3.12b) 


A1  =  A  sin(pnc/a)  (3.12c) 

Equations  (3.8-11)  may  be  cast  into  more  easily  recognizable 
form  by  adding  >z?A  Subtracting  multiples  of  the  equations.  The 
resultant  equations  are 

ao 

Z  A»(— i— )  -  A  *  ( — -— )  +  (-l)m+1  c2yC>  =  0  (3.13) 


,  n'a  +y 
n=l  n  ‘m 


a  -y 
p  ro 


m  m 


i(r-V+  r~“)  -  A'( 


E  A  .  , 

,  n  a  -y  a  +y 
n=l  n  'm  n  'm 


1  +  — ^— )  =  0 


a  +y  a  -y 

p  m  p  m 


(3.14) 


:  -  A'vV*  -  *  0 

p  m 


.  ‘n  a  -i-3 
n=l  n  m 


(3.15) 
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E 

n=l 


A'(- 
n  a 


m 


-B 

n  in 


o  +8 
n  m 


■)  -  A‘( 


a  +B 

p  m 


m 


a  -B 
d  m 


) 


+  6pb2  B  B/dw  =  0  (3.16) 

m  p 

These  equations  are  valid  for  each  m  =  1,  2,  3,  ....  These 

simultaneous  infinite  sets  of  equations  for  the  (A’},  (B  },  and 

n  m 

{Cm>  are  in  the.  general  form  discussed  in  Section  II,  hence  approximate 
solution  by  the  modified  residue-calculus  technique  is  possible. 

Because  of  the  wide  variety  of  possible  variations  on  the  basic 
geometry,  the  solution  of  the  sets  of  equations  will  be  discussed  in 
general  and  then  restricted  cases  of  greater  interest  will  be  treated 
in  detail. 


Solution  of  the  Sets  of  Equations 

Because  of  the  linearity  of  the  system,  it  suffices  to  treat 

excitation  from  region  (A)  and  region  (B)  separately.  Excitation 

from  region  (A)  is  most  applicable  to  specific  cases  of  interest,  so 

the  discussion  of  the  general  solution  will  be  in  that  context.  Suppose, 

then,  that  B=0  and  hence  B=0  in  (3.16), 

The  inhomogeneous  bifurcation  is  a  generalization  of  the  simple 

bifurcation  in  which  e,  =  e  =1.  In  the  simple  case,  £  =  A  =0, 

be  r  *  m  m  * 

so  that  except  for  the  factor  A',  equations  (3.14)  and  (3.16)  reduce 

to  (2.2)  and  (2.3),  where  the  (A'}  are  identified  with  the  (A  }.  The 

n  n 

solution  is  exactly  as  described  in  Section  II  except  that  the  mero- 

morphic  function  g(w)  is  n<  alized  to  give  R  (-a  )  =  A’.  Solution  of 

g  P 

the  general  case  is  based  on  modifications  of  the  function  g(w)  of 
(2.4). 


I 

1*1 


Suppose  that  a  meromorphic  function  h(w)  can  be  constructed  which 


satisfies : 


(h:l)  h(w)  has  simple  poles  at  w  =  - a and  at  w  -  a^,  a^,  a^,  .. 

(n:2)  h(g  )  +  £  h(-8  )  =  0,  a  =  1,  2,  3*  _ 

m  in  o 

(h:3)  h(y  )  +  \  h(-y  )  =  0,  n  =  1,  2,  3,  .... 
m  m  m 

(h:4)  h(w)  ^  Kw~V,  4/3  £  v  <_  5/3,  as  Jwj  -*■  <*•.  (The  exact  value 
of  v  is  determined  by  the  edge  condition  and  will  be 
discussed  subsequently.) 

Consider  the  following  integrals  in  the  complex  w-plane,  where  C  is  a 
closed  contour  enclosing  all  of  the  poles  and  zeros  of  h(w). 


t  CO 

4  (£M) dw  =  £  R.  (a  )(— i~)  -  Rj-a  )(— ^— ) 

7  w+y  ,  n  n  a  +y  n  p  a  -y 

i  ‘m  n=l  n  *m  r  d  'm 


2*j  y  «nB 


t  h(**y  )  =  0j  in  -  1)  2)  3 }  ••••  (3.17) 

m 


1  I  (h(w)  A  h(w) 

27rj  J  w-Ym  n  W+Y_ 


w-Y  m  w+y 

m  'm  n=l 


00  ^ 

)dw  =  E  R.  (a  )(— L_  +  _E__) 
,  n  n  a  -y  a  +y 
n-  'm  n  m 


-R.  (-a  )(•— —  +  — — )  +  h(y  )  +  A  h(-y  )  =  0, 
n  p  a  +y  a  -y  ‘m  m  m  * 

r  p  m  p  m 

m  =  1,  2,  3,  ...  (3.18) 


■k  \  5  W'rV  -  Rh(-“P)(rV) 

4  m  n=i  r>  m  *  n  m 


+  h(-0  )  =  0,  m  =  1,  2.  3,  ...  (3.19) 
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1  i  /b(w)  j.  r  _  Z  p  (n  ir  1  Sn  i 

J  (wIBl  Cp  w+Sn)dW  "  ?.  V“n)(a  -3  *  a  +3  } 

J  £  m  m  n=l  n  m  n  m 


-  R.  (-a  )(— +  — V")  +  M&J  +  5nh(-Bm)  =  °» 
n  Da  +3  a  —3  m  mm 


DP  DP 


p  =  1,  2,  3,  ...  (3.20) 

When  h(w)  is  normalized  to  make  R^C-a^)  =  A',  (3.17-20)  may  be 
identified  term  by  term  with  (3.13-16)  to  generate  the  solutions  to 
the  sets  of  equations.  The  results  are: 


V  -  yan),  n  =  1,  2,  3, 


(3.21) 


C  =  (-l)m+1  -f5-h(-Y _),  m  =  1,  2,  3,  .... 
m  c2  p 

'm 


(3.22) 


B  =  -  -~-h(-3  ),  p  =  1,  2,  3, 
m  b2g  m 

m 


(3.23) 


The  required  function  h(w)  takes  the  form 


»<”>  ■ *  «»«•> 


(3.24) 


The  zeros  { 3 ' }  and  (y * }  reduce  to  the  zeros  (3  )  and  {y  }  of  g(w)  when 
n  n  n  n 

Xm  =  =  0  (equivalently,  =  ec  =  1).  The  constant  L  is  given  in 

(2.6)  and  the  factor  K  is  chosen  to  satisfy  the  normalization  require¬ 
ment.  The  nominal  algebraic  behavior  of  h(w)  is  asymptotically 
-3/2 

w  ,  but  corrections  are  made  by  the  asymptotic  shift  of  3^  and  y^ 
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from  their  nominal  values  of  and  y  .  Properties  of  the  sets  of 
zeros  are  determined  concurrently  by  (h:2),  (h:3) ,  and  the  edge 
condition. 

The  solution  for  excitation  from  region  (B)  proceeds  almost 
exactly  as  for  excitation  from  (A).  The  function  h(w)  differs  in 
that  the  pole  at  -o^  and  the  zero  at  g^  are  eliminated,  so  that  h(w) 
maintains  the  same  algebraic  behavior  as  before.  Subsequent  arguments 
pertaining  to  the  edge  condition  and  asymptotic  behavior  carry  over 
verbatim. 


Edge  Condition  and  Asymptotic  Behavior  of  Zeros 

So  long  as  simple  dielectrics  (e^  <  <»,  ec  <  »)  fill  regions  (B) 
and  (C)  and  only  TE  excitation  is  considered,  the  inhomogeneously- 


filled  waveguide  bifurcation  has  the  same  edge  behavior  as  the  simple 

-3/2 

bifurcation.  Thus,  h(w)  'v  Kw  for  both  cases.  However,  two 
infinite  sets  of  zeros  must  be  determined  simultaneously  to  satisfy 
(h:2)  and  (h:3).  As  in  Section  II,  it  is  possible  to  extract 
information  about  the  asymptotic  zero  positions  by  studying  conditions 
(h:2)  and  (h:3)  for  large  indices. 

Suppose  that  for  n  very  large,  the  propagation  constants  may  be 


written 


“n  =  -jnjr/a 


g  =  -jntt/b 
n 


Yn  =  -jntr/c 


g^  =  -jmr/b  +  S 


35 


Y^  =  -jnit/c  +  D 


where  S  and  D  are  the  assumed  asymptotic  shifts  of  8^  and  y^  from  the 
unshifted  positions.  Imposing  (h:2)  for  m  large  in  the  form 


h(B  )/h(-8  )  =  -  5  yields 
m  m  ro 


jmTi/b+a 


nd-  -Tfe’A) 

_ -jrai/btS 


( — r — /r,*D-)  exD(-2jmnL/b)  - ■■■ - 

-}nnr/b+a  *  J  timn/b  . 

p  V  -jn7i/b+S; 


nu-  -ds^UiKi-  i=£) 

_ -]ror/c+D _ -jmr/a 

aci-  ti™/*  ad-  4=^, 

-jnn/c+D  -jrot/a 


(3.25) 


The  exponential  convergence  factors  in  the  products  have  been  omitted 
to  simplify  the  appearance  of  the  equation.  When  the  products  are 
replaced  by  their  simplified  forms  from  (2.20a)  and  (2.20b)  and  when 
the  quotient  containing  is  replaced  by  its  limiting  value,  equation 
(3.25)  reduces  to 


sin(miT-jSb  )sin(mnc/b-jDc) 
sin(mita/b) 


=  |  C 
2  m 


(3.26) 


Since  a  =  b  +  c,  this  expression  may  be  further  reduced  to 


sin( jSb)sin(mnc/b-jDc) 
sin(imrc/b) 


(3  27) 


The  corresponding  result  for  (h:3)  may  be  obtained  from  (3.27)  upon 
interchanging  b  and  c,  S  and  D,  and  replacing  5^  by  Xm< 
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s in ( j Dc ) s in (mirb/cXj Sb )  1  ,  /,»  0 a 

~  sin(mub/c)  2  Am 

Tor  simple  dielectrics,  both  5  -*■  0  and  X  +  0  as  m  +  ®.  Thus, 

'mm 

the  R.H.S.  of  (3.27)  and  (3.28)  tend  to  zero,  and  hence  so  must  the 

L.H.S.  if  the  equalities  are  to  be  maintained  in  the  limit.  This 

requirement  is  satisfied  by  D  =  S  =  0,  so  that  asymptotically  the 

new  zeros  0’  and  y*  tend  to  the  old  zeros  $  and  y  .  The  choice 
n  n  n  n 

S  =  D.  =  0  is  consistent  with  the  edge  condition,  for  no  algebraic 
correction  to  the  asymptotic  behavior  of  h(w)  is  permitted. 

In  the  most  general  case,  i.e.,  when  one  or  both  dielectrics 
assume  limiting  values  as  electric  or  magnetic  conductors,  such 
direct  limiting  arguments  for  (3.27)  and  (3.28)  fail.  They  fail 
because  the  limits  as  m  -*■  °°  of  the  sine  ratios  fail  to  exist  even 
though  the  R.H.S.  cf  each  equation  does  have  a  limiting  value.  Two 
cases  of  particular  interest,  namely  the  waveguide  step  and  the 
diaphragm,  may  be  solved  as  special  limiting  cases  of  the  inhomogen- 
eously- filled  bifurcation.  In  those  cases,  limit  difficulties  with 
respect  to  asymptotic  behavior  of  the  zeros  reflected  in  (3.27)  and 
(3.28)  can  be  circumvented  by  direct  imposition  of  the  edge  condition. 
Although  the  two  problems  will  be  discussed  in  greater  detail 
subsequently,  the  behavior  of  the  set  of  zeros  of  h(w)  for  their 
solution  will  be  treated  as  special  cases. 

Case  1.  When  the  dielectric  in  region  (C)  is  allowed  to  become 
a  perfect  electric  conductor  and  the  dielectric  in  (B)  remains  simple, 
the  geometry  reduces  to  the  waveguide  step  discontinuity  of  Figure  3.3. 


In  this  case 


,  5  -*■  0  as  a  -*■  «,  but  X  =  -  1  for  each  n.  The  value  of 

D  *3 

X  is  readily  interpreted  as  the  n'"n  rode  reflection  coefficient  at  the 
□ 

interface  between  regions  (C*)  and  (C)  in  the  auxiliary  georetry  of 
Figure  3.2.  Naturally,  when  the  dielectric  be cores  a  perfect  electric 
conductor,  each  node  reflection  coefficient  becores  -1. 

The  limit  of  the  R.H.S.  of  (3.27)  is  zero,  but  the  R.H.S.  of 
(3.23)  is  +  1/2  for  all  a.  The  choice  S  =  0  reduces  the  L.H.S.  of 
(3.27)  to  zero,  so  that  (3.27)  is  satisfied  in  the  liait.  However,  the 
choice  S  =  0  reduces  (3.28)  siaply  to 

sin(jDc)  =  ^  (3.29) 

which  then  implies  that 

2  =  -  %  (3.30) 

Thus  the  zeros  {y^}  are  asymptotically  shifted  from  {y^}  by  the 

amount  D  =  -  but  the  asymptotic  shift  of  the  {S'}  is  zero. 

The  edge  condition  for  the  waveguide  step  is  that  of  a  right- 

angle  perfectly  conducting  wedge.  The  field  behavior  requires  that 
-5/3 

h(w)  'v*  Kw  ,  which  differs  from  the  behavior  for  simple  dielectrics 
by  the  factor  of  w  ^6.  It  would  appear  at  this  point  that  a  new 
form  for  h(w)  must  be  found  to  account  for  the  new  edge  condition. 
However,  the  asymptotic  shift  of  the  zeros  {y^}  automatically 
incorporates  the  algebraic  correction  to  h(w). 

The  ratio  r(w)  of  h(w)  with  shifted  zeros  to  h(w)  with  unshifted 


zeros  is  essentially 


JI(l-w/(y_+D)) 
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However,  substitution  of  asymptotic  values  for  yn  and  reference  to 
(2.20a)  and  (2.20b)  yield  a  simpler  expression  for  the  ratio,  namely 


(  x  _  r(n- jpc/a  )r(i-jwc/a ) 
rdrjDc/a-jwc/s  )F(1) 


(3.32) 


Prom  the  asymptotic  value  of  the  ratio  of  gamma  functions  (T(x+b)/ 
F(x)  i  x*1  as  x  ’  «)  given  in  Noble  [1358],  the  ratio  (3.32)  becomes 


?(w)  ^  K(-j”) 


.Dc 

-3— 


(3.33) 


Substitution  of  the  asymptotic  value  of  D  =  -  jir/6c  yields 

r(w)  '«<  K'w"1/6  (3.34) 

Thus  the  proper  algebraic  correction  to  h(w)  is  automatically  provided. 

Case  2.  When  region  (C)  is  filled  with  a  perfect  electric 

conductor  and  region  (E)  becomes  a  perfect  magnetic  conductor,  the 

geometry  reduces  to  a  waveguide  terminated  by  a  combined  magnetic 

and  electric  wall,  as  in  Figure  3.4.  Although  this  problem  is 

essentially  non-physical  because  magnetic  conductors  fail  to  exist, 

it  is  indeed  useful  as  an  auxiliary  geometry  for  solving  the  problem 

of  an  inductive  diaphragm  in  a  waveguide.  Its  application  to  the 

diaphragm  will  be  discussed  subsequently. 

The  magnetic  conductor  in  region  (B)  gives  the  value  +1  for  each 

£  .  As  in  Case  1,  X  =  -i  for  each  m.  Unlike  Case  1,  however,  the 
mm 

limiting  procedure  applied  to  (3.27)  and  (3.28)  fails  completely, 
indicating  that  the  zeros  fail  to  have  a  specific  asymptotic  shift. 


except  possibly  in  an  average  sense.  Fortunately,  though,  the  edge 

condition  for  Case  2  is  known  exactly,  so  that  edge  condition 

information  can  be  used  to  advantage. 

Because  the  geometry  of  Case  2  may  be  directly  related  to  the 

inductive  diaphragm,  the  edge  condition  is  that  of  the  diaphragm, 

-3/2 

which  like  the  simple  bifurcation,  requires  h(w)  'v  Kw  .  Since 
no  algebraic  correction  to  n(w)  is  required,  the  edge  condition  alone 
yields 


S  =  0 

(3.35a) 

D  =  0 

(3.35b) 

for  the  asymptotic  shifts  of  the  zeros.  Since  the  conditions  (h:2) 
and  (h:3)  for  the  magnetic-electric  wall  case  are  somewhat  stronger, 
it  may  very  well  require  longer  for  the  actual  zeros  to  approach 
the  asymptotic  values.  However,  because  the  actual  shifts  of  the 
first  few  zeros  of  h(w)  are  found  explicitly,  the  numerical  results 
for  scattering  matrix  elements  for  low-order  modes  are  quite 
insensitive  to  ex-rors  in  the  location  of  large-order  zeros  of  h(w), 
hence  in  S  and  D. 

The  essential  test  of  the  proper  c  'oice  of  asymptotic  shifts  is 
the  satisfaction  of  (h:2)  and  (h:3).  Of  course,  when  the  conditions 
assume  proper  limiting  forms,  the  problem  is  greatly  simplified. 
However,  arguments  such  as  those  for  Case  2  may  be  invoked  when  direct 
limiting  procedures  fail.  In  some  cases,  moreover,  explicit  computa¬ 
tion  of  exact  zero  locations  may  lead  to  a  numerical  determination 
of  the  asymptotic  shift  of  large-index  zeros.  In  all,  any  reasonable 


means  that  yields  satisfaction  of  (h:2)  and  (h:3)  and  the  edge 
condition  may  be  useful. 

Construction  of  h(w) 

Since  the  functional  form  of  h(w)  is  known  exactly,  construction 

of' h(w)  consists  in  determining  the  unknown  zeros.  As  in  the  problem 

discussed  in  Section  II, "it  is  not  possible  to  determine  all  of  the 

zeros  exactly.  Again  it  suffices  to  substitute  asymptotic  values  for 

all  but  the  first  few  values  of  0'  and  y'  and  to  find  the  first  values 

n  n 

explicitly.  Suppose  that  M_  values  of  0'  and  M-,  values  of  y'  are  to 

d  n  c*  n 

be  found  explicitly,  and  that  M  =  M_  +  M_.  Actually,  a  polynomial  of 

D  L 

order  M  is  determined  since  exact  location  of  the  zeros  is  errelevant. 
The  function  may  be  written 


h(w)  =  K  P(w }  exp(Lw)n^ft^(w,0' )n^MQ^(w,yt ) 


(w+ap)II(w,a) 


(3.36) 


where  Jr*  B‘'(w,0' )  denotes  omission  of  the  first  M,.  factors  and 

D 

replacement  of  p'  by  0  S  for  n  >  Mn,  and  similarly  for  n^C^(w5y). 


n  "  n 


P(w)  is  a  polynomial  of  order  M.  For  simplicity,  write 


h(w)  =  K  P(w)h'(w) 


(3.37) 


where  h; (v)  incorporates  all  of  the  known  parts  of  h(w). 

Although  the  method  of  determining  directly  the  polynomial 
coefficients  which  is  described  in  Section  II  is  perfectly  applicable 
to  this  case,  an  alternative  method  of  specifying  Hw)  is  superior 
for  numerical  work,  albeit  .tore  complicated.  The  method  consists  of 
writing  P(w)  as  a  known  polynomial  of  order  M  plus  a  correction 


I 


polynomial  of  order  M-l.  For  convenience  in  notation,  let 


Wi  =  1  =  ly  2*  *•**  HB 


M.  =  Yij  i  =  Mg+1,  Hg+2,  ...»  Hs  +  «c  =  M 
B 

=  5^ ,  i  =  1,  2,  Mg 


ti+M  =  V  1  =  mb  +  hb  +  2>  •••»  H 

3 


Thus,  conditions  (h:2)  and  (h:3)  for  the  explicitly  shifted  zeros  read: 


MWJ  +  t.hC-VL)  =  0,  i  =  1,  2,  ....  H 


(3.38) 


Now  write  P(w)  in  the  following  form: 


M  MM  w-W. 

n  (l-w/w..)  +  z  Cf.  n  i-)] 

•  1  i  i  •  W*“n* 

1=1  1=1  3=1  l  3 


(3.39) 


The  first  term  of  P(w)  is  readily  recognizable  as  the  product  of 
unshifted  zeros.  The  second  term,  in  which  the  coefficients  F^, 
i  =  1,  2,  ...,  M,  are  to  be  found,  forms  the  correction  polynomial. 

The  normalization  of  the  second  term  is  such  that  P(W.)  =  F.,  i  =  1,  2, 

li 


...,  M. 


Substituting  the  new  form  of  P(w)  into  (3.38)  yields 


K  n ' (w  )F  t  t  Kh'(H  )  {  II  (1+W  /W.) 
mm  m  m  .  ,  mi 

i=J 


f 
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K  M  -K  -W. 

E  CF  n  (rr^-)]} 

1=1  3=1  i  3 


=  0,  ni  =  lv  2,  3,  . ..,  M 


(3.40) 


Since  the  only  unknowns  are  the  {F^},  which  appear  linearly,  (3.40) 
is  actually  a  set  of  M  simultaneous  linear  equations  for  the  {F^}. 
Solution  of  such  sets  by  numerical  methods  is  well-known.  In  actual 
practice.  Gauss- Jordan  elimination  in  double-precision  complex  arithmetic 
was  used  on  an  IBM-7094  digital  computer. 

The  advantage  of  this  scheme  for  finding  the  polynomial  is  that 
the  unknowns  are  of  first-order,  while  the  unknown  polynomial 
coefficients  in  the  direct  method  are  of  zero-order.  The  alternative 
method  reduces  to  F^  =  0,  i  =  1,  2,  3,  ...,  M,  when  the  geometry 
reduces  to  a  simple  bifurcation;  for  the  same  case,  the  direct  method 
yields  rather  large  values  for  the  unknowns.  The  best  test  of  the 
polynomial  determination  is  the  satisfaction  of  (h:2)  and  (h:3),  or 
equivalently,  (3.38).  The  correction  polynomial  method  is  consistently 
superior  to  the  direct  polynomial  method  in  satisfying  (3.38) 
numerically . 


Application  to  Specific  Examples 

Simple  dielectric  fillings.  Apart  from  the  homogeneous 
bifurcation,  the  simplest  geometry  is  that  of  Figure  3.1,  in  which 
regions  (B)  and  (C)  are  filled  with  simple  dielectrics.  For  a  fixed 
waveguide  width  a,  the  reflection  coefficient  varies  with  the  position 
of  the  septum  (i.e.,  with  the  ratio  c/a).  Over  intervals  in  which 
both  guide  (B)  and  guide  (C)  are  cut  off,  the  magnitude  of  the 


reflection  coefficient  is  unity.  Because  the  cut-off  dimension  also 
varies  with  the  nature  of  the  dielectric  filling,  the  interval  of 
total  reflection  varies  with  the  dielectric  constant.  Figure  3.5 
illustrates  the  cut  off  phenomenon  for  a  singly-filled  bifurcation. 
Region  (C)  is  filled  with  a  dielectric  ec,  but  region  (B)  is  unfilled. 
For  the  guide  dimension  a  =  0.75A,  the  reflection  coefficient  for 
the  TE,q  mode  from  (A)  is  given  for  £c  =  1.0,  2.0,  and  3.0.  Note 
the  coincidence  of  the  reflection  coefficients  up  to  the  point 
c/a  -  0.33  at  which  point  the  unfilled  guide  (B)  cuts  off.  Also  of 
interest  are  the  limiting  values  as  c/a  =  1,  in  which  cases  the 
junction  becomes  that  between  an  unfilled  guide  and  a  filled  guide 
of  identical  dimension.  The  cut  off  points  for  the  guide  (O'  are 
readily  observable  as  break-points  in  the  curves  of  Figure  3.5.  Both 
magnitude  and  phase  of  the  reflection  coefficient  for  discrete  values 
of  c/a  for  the  same  geometry  are  giver  in  Table  3.1. 

Reflection  coefficients  for  the  doubly-filled  bifurcation  are 
depicted  in  Figure  3.6.  In  each  case,  guide  (B)  is  filled  with  a 
dielectric  =  3.0,  while  results  are  given  for  ec  =  1.0,  2.0,  and 
3.0.  Again  break-points  in  the  curves  may  be  identified  with  cut  off 
dimensions  in  the  waveguide.  As  before,  a  =  0.75A. 

Step  discontinuity.  Of  greater  practical  importance  than  the 
simple  dielectric-filled  bifurcation  is  the  waveguide  step  discontinuity 
of  Figure  3.3.  Extensive  results  for  .step  discontinuities  of  various 
kinds  are  given  in  Marcuvitz  [1951],  where  primarily  variational 
techniques  are  employed  to  finu  the  dominant  mode  reflection 
coefficient  and  equivalent  circuit  reactance.  However,  no  results 


TABLE  3.1 


Reflection  Coefficient  for  Inhomogeneously-Filled  Bifurcated  Waveguide 


Reflecticn 

Coefficient  (Magnitude, 

Angle) 

e  =1.0,  e  =1.0 

e  =2.0,  e  =1.0 

e  =3.0,  e.  =1.0 

0.0 

0. 

0. 

0. 

0.1 

(0.051,  36.3°) 

(0.051,  36.1°) 

(0.051,  36.1°) 

0.2 

(0.148,  50.6°) 

(0.148,  50.5°) 

(0.147,  50.3°) 

0.3 

(0.420,  53.6°) 

(0.420,  53.1°) 

(0.416,  32.5°) 

0.4 

(1.000,  112.3°) 

(1.000,  109.1°) 

(0.847,  100.6°) 

0.5 

(1.000,  127.2°) 

(0.656,  110.8°) 

(0.620,  134.3°) 

0.6 

(1.000,  112.3°) 

(0.417,  135.5°) 

(0.510,  152.0°) 

0.7 

(0.420,  53.6°) 

(0.318,  151.8°) 

(0.433,  162.5°) 

0.8 

(0.148,  50.6°) 

(0.264,  164.5°) 

(0.371,  170.5°) 

0.9 

(0.051,  36.3°) 

(0.244,  174.8°) 

(0.358,  177.7°) 

1.0 

0. 

(0.252,  180.0°) 

(0.365,  180.0°) 

seem  to  be  given  for  the  junction  of  guides  of  different  dimension 
and  different  dielectric  filling ,  which  Is  the  most  general  case 
of  Figure  3.3* 

Because  the  step  discontinuity  is  merely  an  extension  of  the 

inhomogeneously-filled  bifurcation,  scattering  matrix  representations 

of  the  discontinuity  can  be  readily  found  by  the  MRCT.  Although  it 

often  suffices  to  determine  only  the  dominant  mode  reflection 

coefficients,  knowledge  of  the  scatcering  matrices,  even  though 

truncated  at  a  small  size  (C'S,  for  example),  makes  possible  very 

simple  solutions  of  each  problems  as  the  thick-wall  phased  array 

(Section  IV).  In  addition,  knowledge  of  the  high-order  mode 

coefficients  for  the  scattered  field  permit  very  accurate  determination 

of  the  electric  or  magnetic  field  at  the  discontinuity.  An  example 

AA 

of  the  truncated  scattering  matrix  S  for  a  step  discontinuity  is 
given  In  Table  3.2.  The  guide  width  a  =  0.75X  and  c/a  =  0.3;  region 
(B)  is  unfilled. 

Even  though  only  «.  very  small  numoer  of  zeros  of  the  function 
h(w)  are  determined  explicitly  in  the  solution  of  the  step  problem, 
availability  of  approximate  values  for  the  high-order  mode  coefficients 
permits  computation  of  the  electric  field  at  the  plane  z  =  0.  Assume 
that  a  TE^q  mode  of  unit  amplitude  is  incident  from  region  (A).  The 
coefficients  of  the  first  twenty  reflected  modes  are  found  by  the 
tiT.CT  and  the  total  field  is  calculated  from  the  truncated  Fourier 
series  at  z  =  0.  The  amplitudes  of  total  electric  field  for  step 
widths  of  0.2a,  0.3a,  0.4a,  and  0.5a  are  shown  in  Figure  3.7.  It  is 
important  to  note  that  the  fields  are  computed  in  region  (A),  so  that 


Truncated  Scattering  Matrix  S  for  Waveguide  Step  Discontinuity 
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the  extremely  small  value  (almost  zero)  of  the  field  free  x  =  0  to 
x  =  c  for  each  case  illustrates  the  satisfaction  of  the  requirement. 

=  0onx  =  0tox  =  c  for  even  the  truncated  series.  Comparison 
of  the  fields  in  the  aperture  c  <  x<aas  computed  in  region  (A)  with 
the  corresponding  values  computed  in  region  .(B)  have  shewn  agreement 
to  several  significant  figures.  Although  the  mode  coefficients  found 
by  the  KRCT  are  only  approximate,  the  field  found  froa  theft  is  probably 
very  near  to  a  least-neai. -square  approximation  to  the  actual  value 
because  of  the  Fourier-Bessel  inequality. 

The  reflection  coefficient  as  a  function  of  relative  step  width 
c/a  is  shown  in  Figure  3.8.  Region  (B)  is  filled  with  dielectric 
=  1.0,.  2.0,  and  3.0.  The  guide  width  a  =  0.75X.  The  limiting 
case  c/a  =  0  yields  the  reflection  coefficient  for  the  air-dielectric 
interface  in  a  waveguide;  the  case  c/a  =  1  yields  the  reflection 
coefficient  -1  for. a  shorted  waveguide.  As  previously,  the  break 
points  in  the  magnitude  and  phase  curves  occur  as  guide  (B)  cuts  off. 

Magnetic  and  electric  wall.  The  magnetic  wall-electric  wall 
termination  of  a  waveguide  illustrated  in  Figure  3.4  may  be  solved 
by  the  MRCT  as  discussed  under  Case  2.  The  geometry  of  Figure  3.4 
is  useful  as  an  auxiliary  geometry  in  the  solution  of  the  waveguide 
diaphragm  of  Figure  3.9.  An  additional  auxiliary  geometry  is  the 
shorted  waveguide,  the  solution  of  which  is  trivial. 

Considei  equal  excitation  oF  the  diaphragm  of  Figure  3.9  from 
regions  (A)  and  (A1)  in  a  symmetric  sense,  i.e.,  with  electric  vectors 
in  both  regions  pointing  in  the  same  direction.  Then  at  the  plane 
z  =  Q,  the  electric  vectors  add  in  the  aperture  and  vanish  on  the 
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diaphragm,  but  the  transverse  magnetic  vectors  cancel  exactly  in  the 
aperture.  Thus  placement  of  a  magnetic  wall  in  the  aperture  has  no 
effect  on  the  fields,  and  so  far  as  region  (A)  is  concerned,  the 
geometry  is  that  of  Figure  3.4.  Alternative  equal  but  antisymmetric 
excitation  of  the  diaphragm  from  (A)  and  (A’}  causes  the  total  electric 
vector  to  vanish  on  the  plane  z  =  0,  so  that  an  electric  wall  may  be 
placed  there  without  effect.  Then,  the  geometry  for  region.  (A)  reduces 
to  a  shorted  waveguide. 

Suppose,  now,  that  both  symmetric  and  antisymmetric  excitation 

are  applied  to  the  diaphragm  simultaneously,  so  that  addition  occurs 

AA 

in  (A)  but  exact  cancellation  occurs  in  (A').  Let  be  the 

scattering  matrix  determined  approximately  for  the  magnetic-electric 

wall,  and  let  sf  =  -  I  be  the  scattering  matrix  for  the  shorted 

AA 

waveguide.  Then  the  resultant  scattering  matrix  S  for  the  diaphragm 
is  given  by 

gAA  =  l^AA  +  SAA}  =  (3.41) 

If  R^  is  the  complex  reflection  coefficient  for  the  dominant  mode  for 
the  magnetic-electric  wall,  then  the  reflection  coefficient.  R  for 
the  diap’.ii-agm  of  the  same  dimension  is 

R  =  ■|(R1-1)  (3.42) 

Because  no  power  is  transmitted  beyond  the  magnetic-electric  wall,  the 
magnitude  of  is  always  exactly  unity.  The  phase  changes  from  0  to 
tt  as  the  ratio  c/a  goes  from  0  to  1.  The  magnitude  and  phase  of  R  as 
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a  function  of  c/a  for  the  diaphragm  is  given  in  Figure  3.10.  The  guide 
width  is  a  =  0.75A. 

An  additional  application  of  the  scattering  matrix  description  of 
the  magnetic-electric  wall  geometry  is  diffraction  by  an  infinite  array 
of  thin  conducting  strips.  This  problem  is  discussed  in  Section  V. 

Alternative  Solution  Methods 

By  the  very  nature  of  the  geometry  of  the  problems  related  to  the 
simple  waveguide  bifurcation,  it  is  apparent  that  alternative  methods 
of  solution  exist.  The  generalized  scattering  matrix  procedure  is 
directly  applicable  fo  many  of  the  problems,  since  the  geometries  can 
often  be  reduced  to  combinations  of  exactly  solvable  junctions. 
Variational  and  integral  equation  methods  such  as  those  described  by 
Marcuvitz  [1951]  have  been  applied  with  success  to  many  such  problems. 
In  any  specific  instance,  the  choice  of  method  of  solution  can  be  made 
in  many  ways. 

As  applied  to  geometries  related  to  the  waveguide  bifurcation, 
the  modified  residue-calculus  technique  possesses  some  great  advantages. 
First,  it  is  extremely  general.  All  of  the  problems  discussed  in 
this  section  were  solved  by  use  of  the  same  computer  program,  whether 
the  problem  was  the  simple  bifurcation  to  be  solved  exactly  or  the 
magnetic-electric  wall.  Second,  it  generates  scattering  matrix 
representations  of  discontinuities,  such  as  those  obtained  by  the 
multiple-reflection,  Neumann-series  method.  On  the  other  hand,  it 
can  generate  the  first  M  scattered  mode  coefficients  for  the  first 
incident  mode  without  resorting  to  an  matrix  inversion.  Inversion 


58 


of  large-order  matrices  is  a  major  drawback  to  some  methods  in  the 
solution  of  integral  equations.  Not  least  important  is  the  automatic 
satisfaction  of  the  edge  condition.  In  both  direct  matrix  inversion 
methods:  as  well  as  the  generalized  scattering  matrix  method,  it  is 
difficult  to  prove  directly  that  the  edge  conditio, a  is  satisfied. 

In  the  MRCT,  satisfaction  is  a  priori  guaranteed. 
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IV.  THICK-WALL  PHASED  ARRAY 

At  the  present  time,  there  is  great  interest  in  antennas  which 
arc  physically  stationary  but  which  electronically  scan  a  relatively 
narrow  beam  throughout  a  region  of  space.  Two  particular  applications 
are  radio  telescopes,  in  which  the  antenna  is  often  so  large  that 
mechanical  steering  is  prohibitively  difficult,  and  tracking  radars, 
in  which  a'  physical  antenna  has  too  great  inertia  to  be  rapidly  and 
repeatedly  scanned  over  a  large  area  while  tracking  several  fast- 
moving  targets.  In  electronically  scanned  antennas,  phase  changes 
between  elements  steer  the  antenna  beam  without  need  for  mechanical 
positioning . 

One  type  of  electronically  scanned  antenna  is  a  phased  array  of 
open-end  rectangular  waveguides,  a  segment  of  which  is  shown  in 
Figure  4.1.  Although  the  number  of  elements  (waveguides)  in  such  an 
array  is  necessarily  finite,  it  is  not  uncommon  in  practice  for  the 
number  of  elements  to  be  extremely  large.  For  purposes  of  analysis,  it 
is  convenient  to  assume  that  the  array  is  infinite  in  extent  and 
periodic.  The  infinite  array  approximation  is  good  for  the  majority 
of  elements  away  from  the  edges  of  the  finite  array.  It  is  further 
assumed  that  the  waveguides  are  excited  by  the  dominant  TE  mode  with 
uniform  amplitude  and  periodic  phasing  to  scan  the  beam.  To  be 
determined  for  the  array  are  the  aperture  fields,  from  which  the 
radiation  pattern  may  be  computed,  and  the  reflection  coefficient  for 
thi'  dominant  mode  in  each  waveguide  as  a  function  of  scan  angle. 

The  infinite  rectangular  array  was  discussed  by  Wu  and  Galindo 
[1966]  for  H-plane  and  quasi-E-pIane  scan  for  infinitely  thin 


Figure  4.1  Phased  array  of  open-end  rectangular  waveguides 
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waveguide  walls  and  by  Galindo  and  Hu  [1965]  for  walls  of  finite 
thickness.  The  case  of  thin  walls  and  H-plane  scan  is  equivalent  to 
diffraction  by  an  infinite  set  of  thin  parallel  plates,  which  nay 
be  solved  exactly  by  the  Viener-Bopf  technique  (Carlson  and  Heins 
[1947 J  or  by  the  conventional  residue-calculus  technique  (Berz  [19513, 
Whitehead  [19513).  However,  the  corresponding  problem  with  walls 
of  non-zero  thickness  cannot  be  solved  exactly.  Radiation  at 
broadside  frca  an  array  of  thick  plates  was  studied  by  Lee  [19573, 
who  presented  a  solution  based  on  the  Wiener-Hopf  technique.  By  a 
numerical  solution  to  an  integral  equation  for  the  aperture  field, 
Galindo  and  Wu  [19553  were  able  to  find  reflection  coefficients  and 
aperture  fields  for  several  cases  of  wall  thickness  and  waveguide 
dissension  for  angles  of  scan  frca  0®  to  90°. 

Not  unexpectedly,  the  results  of  Galindo  and  Wu  indicate  that 
the  reflection  coefficient  solutions  for  the  thin  wall  array  are  not 
very  good  approximations  to  the  reflection  coefficients  for  arrays 
with  walls  of  non-zero  thickness.  Because  a  major  consideration  in 
array  design  is  to  minimize  reflection  at  the  aperture,  it  is 
essential  tc  know  the  variation  of  reflection  coefficient  with  wall 
thickness  and,  in  some  circumstances,  with  the  dielectric  properties 
of  the  material  filling  the  waveguide.  (A  technique  sometimes  used 
to  match  an  array  is  to  fill  the  waveguides  with  dielectric  material.) 
But,  it  is  extremely  useful  to  determine  the  required  information  in 
the  simplest  possible  way,  so  that  as  little  time  as  possible  is 
consumed  in  determining  the  required  array  characteristics. 


Because  of  toe  nature  of  the  array  geometry,  fields  in  each  region 


cay  be  written  in  modal  expansions.  As  a  result,  the  reflection  and 
transmission  characteristics  of  the  array  cay  be  completely  described 
by  a  set  of  scattering  matrices .  Since  the  model  for  the  array  cay 
be  reduced  to  two  auxiliary  junction  problems ,  the  generalized 
scattering  matrix  procedure  may  be  readily  applied  to  the  array 
geccetry.  The  MSCT  solution  of  one  of  the  related  junction  problems 
makes  possible  an  extremely  simple  solution  to  the  array  problem.  In 
particular,  it  permits  determination  of  an  approximate  expression  for 
the  dominant  mode  reflection  coefficient  in  a  fora  suitable  for  simple 
hand  calculation  frcea  data  available  in  tabular  or  graphical  form. 

Thin-wall  Array 

The  solution  for  the  array  with  thin  walls  and  H-plane  scan  forms 
the  basis  for  the  solution  of  the  corresponding  thick-wall  array.  As 
previously  mentioned,  this  problem  was  discussed  in  detail  by  Vu  and 
Galindo  [1956 j,  who  obtained  an  exact  solution  for  dominant  mode 
excitation  by  means  of  the  conventional  residue-calculus  technique. 
Because  solution  in  terms  of  a  generalized  scattering  matrix  is 
essential  to  further  results,  the  solution  for  the  thin-wall  array 
with  general  TE  excitation  will  be  summarized  here. 

The  thin-wall  array  for  H-plane  scan  is  equivalent  to  a  set  of 

thin,  perfectly  conducting  plates  uniformly  spaced  (Figure  4.2). 

Each  parallel  plate  waveguide  thus  formed  is  excited  in  the  region 

z  <  0  by  a  TE  mode  of  unit  amplitude.  The  electric  field  in  the 

p,o 

waveguide  (region  (A))  may  be  written 
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E  =  sin(^-  (x-sa))  exp(-jmu-ja  z) 

y  a  P 


+  E 
n=l 


A  sin{ —  (x-sa))  exo(-jau+ja  z) 
n  a  *  n 


(4.1) 


and  in  the  open  region  (F) 


Z  G  expC-j^x-jr^z) 
n=—» 


(4.2) 


where 


u  =  k  a  sin  8 
o 


£  =  (2nx+u)/a 
n 


a .  =  »4  2-(ns/a)2 
n  o 


Space  modes  for  which  has  a  non-zero  real  part  propagate  away  from 
the  array  and  carry  real  power.  The  number  of  propagating  space  modes 
is  a  function  of  both  guide  width  a  and  scan  angle  9.  Expressions 
proportional  to  the  transverse  magnetic  field  H  are  found  by 
differentiating  (4.1)  and  (4.2)  with  respect  to  z. 

Upon  matching  the  transverse  field  components  at  z  =  0  and 
Fourier  analyzing  the  resultant  expressions  ?>s  described  in  Section 
III,  the  following  infinite  sets  of  equations  are  obtained: 


ft 

Z  n  A  [(-l)n  exp(ju)-l]( — 

.  n  a  -t 

n=l  n  q 

-  p[(-l)P  exp(ju)-l](--^-  )  =  0 

P  q 


(4.3) 
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£  n  A  C(-l)n  exo(ju)-l](-4— ) 

.  n  “  a.  +t 

n=l  n  o 


-  d[C-1)P  exp(ju)-i3  +  2  —  t  G  =  0  (4.4) 

■  o  q 


The  solution  to  (4.3)  and  (4.4)  is  constructed  by  aeans  of  a 
oeroaorphic  function  g(w)  which  takes  the  fora 

g(w)  =  p  exp[j(w+o )  —  £n 


* 


-•  l+o  /a  “>  1-w/t 

x  n  H  -P,-n)  n  (- — ~) 

.  l-w/a  1+a  /r 

n=l  n  n=-»  p  n 


(4.5) 


The  node  coefficients  then  are 


A  = 


R  (a  ) 

-g-JL. 


q  aC(-l)<1  exp(ju)-l] 


'»  ?  ■  3,  ... 


(4.6) 


«  2,2t  ’ 


q  =  0,  +1,  +2,  ... 


(4.7) 


These  solutions  immediately  yield  the  elements  of  the  scattering 
matrices  corresponding  to  excitation  from  the  (A)  or  waveguide  region. 
Since  the  incident  wave  was  the  p  n  mode  of  unit  amplitude,  the 
scattering  matrix  elements  are: 


AA 

S  (q,p)  =  Aq  ,  q  =  1,  2,  3,  ... 

FA 

S  (q,p)  =  Gq  ,  q  =  0,  +  1,  +  2,  . 


(4.8) 

(4.9) 


;  ^ 


where  Aq  and  Gq  are  given  in  (4.6)  and  (4.7),  respectively.  Solution 
for  each  value  of  p  generates  the  entire  scattering  matrices,  which 
are  implicitly  functions  of  the  scan  angle  8. 

The  related  scattering  matrices  corresponding  to  excitation  from 
the  (F)  or  free-space  region  are  constructed  in  a  similar  manner.  Assume 
that  the  p^  free-space  mode  of  unit  amplitude  is  incident  upon  the 
array,  the  electric  fields  in  the  waveguide  region  (A)  may  be  written 
(mth  guide) 

CD 

E  =  Z  A  sin(^kx-ma))exD(-jmu+ja  z)  (4.10) 

v  .  n  a  ~  n 

J  n=l 

and  in  the  open  region  (F) 


Ey  =  exp(-jCpx+jtpz) 


+  Z  G  exp(-]5  x--jT  z) 
n  r  J  n  J  n 

n=-«° 


(4.11) 


Upon  matching  E  and  H  fields,  the  resultant  equations  may  be 
reduced  to 


Z  n  A  C(-l)n  exp(]u)-l]( - )  -  2  —  t  6^=0 

.  n  a-x  upp 

n=l  n  q  v  v 


q  =  0,  +  1,  +  2,  ... 
2 

iT  q  q 


(4.12) 


Z  n  A  C(-lj  exp(ju)-l]( — — — )  +  2  —  t  G  =  0 

n  a+T  "  «  " 

n=l  n  q 


q  =  0,  +  1,  +  2,  ...  (4.13) 
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The  solution  to  this  set  is  constructed  from  a  function 


f(w),  ’’here 


-  2  —  t  e:n>Cj(w-r  )  —  £n  2j 

v  p  -  J  p  it  , 


®  1-7  /a  ®  1-w/t 

*  ”  <i :^rL>  "  ‘rTTT-) 

n=l  n  n=-«  p  n 


(4.14) 


The  mode  coefficients  and  hence  scattering  matrix  elements  are  then 


A  =  SAt(q,p)  =  - -  ,  q  =  1,  2,  3, 

q  q[(-l)q  exp(ju)-l] 


(4.15) 


G  =  SFF(q,p)  = 


7Tf(-T  ) 


»  ^  =  -  >  +.1*  +.2,  ... 


(4.16) 


The  properties  of  the  thin-wall  array  for  each  scan  angle  0  are 

AA  AF  FA 

completely  described  by  the  four  scattering  matrices  S  ,  S  ,  S  , 

FF 

and  S  .  The  dominant  mode  reflection  coefficient,  which  is  equal  to 
AA 

S  (1,1),  as  a  function  of  scan  angle  is  depicted  in  Figure  4.3  for 
three  values  of  waveguide  width  a.  The  break  points  in  the  reflection 
coefficient  curves  occur  when  a  second  space  mode  changes  from 
evanescent  to  propagating. 


Thick-wall  Array 

The  model  for  the  thick-wall  array  for  H-plane  scan  used  by 
Galindo  and  Wu,  and  to  be  used  here,  is  a  set  of  thick,  perfectly 


conducting,  parallel  plates,  as  shown  in  Figure  4.4.  The  spaces, 
between  the  plates  may  be  considered  to  be  parallel  plate  waveguides, 
each  guide  excited  from  the  region  z  <  0  by  a  TE^q  mode.  The 
amplitude  of  excitation  is  uniform  in  each  guide,  and  the  phase 
progresses  from  guide  to  guide  to  scan  the  main  beam  of  the  radiation 
pattern.  The  fields  in  the  open  (F)  region  admit  a  modal  expansion 
of  the  form  (4.2)  because  of  the  periodicity  of  the  structure  and 
the  excitation. 

Despite  the  modal  expansions  for  the  fields  on  either  side  of 
the  interface  z  =  0,  matching  of  the  transverse  components  does  not 
lead  to  a  set  of  equations  of  the  form  solvable  by  the  modified 
residue-calculus  technique.  However,  the  modal  nature  of  tie  fields 
in  regions  (B)  and  (F)  suggests  description  of  the  array  in  terms 
of  generalized  scattering  matrices.  The  scattering  matrices  may 
be  derived  via  the  multiple-reflection,  Neumann-series  method 
described  in  Section  II  for  a  suitable  auxiliary  geometry. 

The  auxiliary  geometry  appropriate  for  the  thick-wall  array  is 

that  of  Figure  4.5.  The  junction  at  z  =  0  is  precisely  that  of 

the  thin-wall  array.  An  exact  scattering  matrix  description  of  that 

junction  was  given  earlier  in  this  section.  Let  the  scattering 

AA  AF 

matrices  for  the  junction  have  the  subscript  1,  i.e.,  , 

FA  FF 

Si  ,  Si  .  On  the  other  hand,  the  junction  at  z=-6  in  each  waveguide 
is  precisely  the  step  discontinuity  described  in  Section  III.  Although 
the  general  waveguide  step  problem  cannot  be  solved  exactly,  a  very 
accurate  truncated  scattering  matrix  description  of  the  junction  can 
be  found  by  the  MRCT ,  as  in  Table  3.2.  Let  the  scattering  matrices 


?2 


for  this  junction  carry  the  subscript  2,  i.e. ,  S"~  ,  ,  S^“,  - 

33  3F 

The  four  scattering  matrices  for  the  composite  array  -  3  ,  S  *, 

T2  ?? 

S*  ,  S~*  -  will  have  no  subscript. 

Following  the  multiple-reflection  phenomena  as  depicted  in 
Figure  2.4  and  reducing  6  to  zero  yields  Neumann  series  representations 
for  the  composite  scattering  matrices.  When  seamed,  these  give 


Str  =  S“  +  S^Cl-sJ*  Sf]"1  S ^  S^?  (4.20) 

It  Bust  be  noted  that  each  scattering  matrix  with  subscript  1  is  an 
implicit  function  of  scan  angle,  so  that  the  composite  scattering 
matrices  also  are  dependent  upon  angle. 

Because  the  matrix  inversion  must  be  performed  numerically  it  is 
important  to  truncate  the  scattering  matrices  at  the  smallest  size 
which  yields  suitably  accurate  results.  Since  the  scattering  matrix 
elements  must  also  be  computed,  keeping  the  truncated  matrices  as  small 
as  possible  also  reduces  computation  time.  Just  what  constitutes 
"suitably  accurate"  is  a  matter  of  judgment  in  each  case,  but  numerical 
results  of  this  section  give  an  indication  of  variation  of  dominant 
mode  reflection  coefficient  for  matrix  sizes  of  1x1,  3x3,  and  5x5. 


One  particular  advantage  of  formulating  the  thick-wall  array  in 
terms  of  the  scattering  matrix  equations  (4*17-20)  is  that  when  the 
wall  thickness  c  reduces  to  zero,  the  resultant  scattering  matrices 


reduce  to  exactly  the  th in-wall  scattering  matrices.  This  result 
occurs  because  for  the  case  c=0,  the  scattering  matrix  description 
of  junction  2  yields 


(4.21) 


(4.22) 


Substitution  of  these  values  into  (4.17-20)  reduce  the  equations  to 
simply  the  corresponding  thin-wall  scattering  matrices. 

Numerical  Results 

The  scattering  matrix  methods  just  described  were  applied  to  a 
variety  of  cases  for  the  thick-wall  array.  Many  of  the  choices  of 
wall  thickness  and  guide  width  were  made  to  coincide  with  those  given 
by  Galindo  and  Hu  C1966]  to  permit  comparison  of  results.  In  each  case, 
dimensions  were  chosen  to  allow  only  one  propagating  mode  in  each 
waveguide. 

In  each  comparison  case,  the  scattering  matrix  method  gave 
numerical  results  for  magnitude  and  phase  of  the  reflection  coefficient 
versus  scan  angle  which  were  indistinguishable  from  the  results  given 
graphically  by  Galindo  and  Hu.  By  way  of  further  comparison,  though, 
the  method  of  Galindo  and  Hu  was  considerably  more  complicated.  For 
each  choice  of  width,  thickness,  and  scan  angle,  the  appropriate 


integral  equation  vas  approximately  solved  by  Gal er kin’s  aetbod,  and 
involved  solution  of  a  30x30  set  of  linear  equations.  Then  a 
variational  correction  vis  applied  to  obtain  the  reflection 
coefficient.  On  the  other  hand,  the  scattering  nstrix  method 
required  solution  of  two  auxiliary  problems  first  by  quite  eleaentary 
means,  and  then  sisple  ccabinaticn  in  the  fora  of  equations  (4.17-20). 
Each  combination  involved  the  numerical  inversion  of  a  5x5  matrix. 
Since  the  solution  to  the  thin-vall  array  problea  is  common  to  all, 
thick-vall  arrays  for  the  sane  width  a,  and  since  the  solution  of 
each  step  width  is  valid  for  all  scan  angles,  the  number  of 
calculations  for  many  related  cases  is  considerably  reduced. 

The  results  of  computations  for  one  choice  of  spacing  (a  =  0.62a) 
and  four  wall  thickness  from  c/a  =  0  to  c/a  =  0.12  are  summarized  in 
Tables  4.1-4.  The  reflection  coefficients  for  the  thin-vall  array 
(Table  4.1)  are  based  upon  the  exact  solution.  In  each  of  the  other 
three  cases,  approximate  results  for  three  sizes  of  truncated 
matrices  are  given.  The  reflection  coefficients  found  by  the  5x5 
scattering  matrix  inversion  agrees  exactly  with  the  values  found  by 
Galindo  and  Hu,  and  will  be  considered  to  be  the  "correct”  value.  The 
results  are  also  depicted  graphically  in  Figure  4.6. 

It  is  significant  to  compare  the  value  obtained  from  the  1x1 
matrix  with  the  correct  value.  Particularly  for  small  scan  angles, 
the  lxl  value  is  very  close  to  the  correct  value,  often  well  within 
one  per  cent  in  magnitude  and  five  degrees  in  phase.  Since  a  lxl 
matrix  is  simply  a  complex  number,  the  computation  of  the  reflection 
coefficient  is  almost  trivial.  It  is  particularly  important  to 
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TABLE  4.1 

Reflection  Coefficient  vs.  Scan  Angle  for  Tfcin-Wall 
Phased  Array  (a/A  =  0.6205) 


Scan  Angle  Reflection  Coefficient  (Magnitude,  Angle) 


0° 

10° 

20° 

30° 

40° 

50° 

60° 

70° 

80c 


(0.2561,  147.5°) 
(0.24S0,  144.5°) 
(0.2268,  133. S°) 
(0.1878,  109.3°) 
(0.0525,  32.0°) 
(0.0021,  32.4°) 
(0.0050,  32.4°) 
(0.0298,  32.1°) 
(0.0745,  31.8°) 
(0.1438,  31.7°) 


90° 
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TABLE  4.  2 

Reflection  Coefficient  vs.  Scan  Angle  for  Thick-Wall  Phased  Array 

(a/X  =  0.6205.  c/a  =  0.02) 


:  i 

Scan  Angle 

! 

i 


Reflection  Coefficient  (Magnitude,  Angle) 


1*1  Matrix  3*3  Matrix  Matrix 


c° 

(0.2730, 

149.7°) 

(0.2740, 

148.6°) 

(0.2743, 

148.0°) 

10° 

(0.2656, 

146.8**) 

(0.2665, 

145.8°) 

(0.2667, 

145.2°) 

20° 

(G.2419, 

137.2°) 

(0.2427, 

136.3°) 

(0.2428, 

135.7°) 

30° 

(0.1974, 

114.9°) 

(0.1978, 

114.0°) 

(0.1977, 

113.4°) 

40° 

(0.0395, 

51.2°) 

(0.0393, 

49.8°) 

(0.0390, 

48.9°) 

50° 

(0.0133, 

167.7°) 

(0.0182, 

168.6°) 

(0.0182, 

159.3°) 

60° 

(0.0164, 

160.6°) 

(0.0162, 

161.6°) 

(0.0162, 

162.4°) 

70° 

(0.0195, 

73.5°) 

(0.0192, 

72.1°) 

(0.0188, 

71.3°) 

30° 

(0.0608, 

44.3°) 

(0.0608, 

43.2°) 

(0.0606, 

42.4°) 

90° 

(0.1295, 

37.8°) 

(0.1298, 

35.8°) 

(0.1297, 

36.1°) 
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TABLE  4.4 

Reflection  Coefficient  vs.  Scan  Angle  for  Thick-Wall  Phased  Array 

(a/A  =  0.6205,  c/a  =  0.12) 


Reflection  Coefficient  (Magnitude,  Angle) 


Scan  Angle  - 

lx?  Matrix 


0° 

(0.4261, 

158.2°) 

10° 

(0.4193, 

156.6°) 

20° 

(0.3961, 

151.5°) 

30° 

(0.3429, 

140.8°) 

40° 

(0.1711, 

146.2°) 

50° 

(0.1902, 

160.5°) 

60° 

(0.18e8, 

159.8°) 

70° 

(0.1783, 

153.0°) 

80° 

(0.1667, 

139.1°) 

90° 

(0.1719, 

115.9°) 

3x3  Matrix  5x5  Matrix 


(0.4259, 

153.3°) 

(0.4251, 

151.4°) 

(0.4175, 

151.9°) 

(0.4157, 

150.0°) 

(0.3894, 

147.2°) 

(0.3842, 

145.5°) 

(0.3253, 

137.3°) 

(0.3140, 

135.9°) 

(0.1553, 

153.6°) 

(0.1494, 

158.2°) 

(0.1898, 

166.7°) 

(0.1898, 

169.6°) 

(0.1876, 

166.1°) 

(0.1873, 

169.1°) 

(0.1697, 

160.1°) 

(0.1668, 

163.9°) 

(0.1441, 

146.0°) 

(0.1349, 

150.9°) 

(0.1306, 

116.3°) 

(0.1097, 

118.7°) 
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notice  that  simple  exact  reflection  and  transmission  coefficients  for 
the  thin-wall  array  cay  be  combined  with  reflection  and  transmission 
coefficients  for  a  waveguide  step  discontinuity  (found  by  the  MRCT 
or  even  from  graphical  results,  e.g.  Marcuvitz  [1951])  to  yield  an 
excellent  approximate  result  for  the  thick-wall  array.  The  simplified 
expression  for  the  array  reflection  coefficient  may  be  written 


BA  AA  AB 

p  _  pBB  A  *2  *1  L2 

2  1  -  r^a  raa 

4  1 


(4.23) 


where  R™  *  sAA(l,l) 

R^  =  S^(l,l) 

Rf  =  S2B(1>1) 

tab  =  SAB(1,1) 

TBA  =  SBA(1,1) 

Although  the  array  of  Figure  4.4  is  shown  without  dielectric 
filling  in  region  (B),  the  solution  for  dielectric  filling  proceeds 
exactly  as  previously  described.  The  corresponding  dielectric-filled 
step  discontinuity  is  used  as  part  of  the  auxiliary  geometry. 

The  mode  coefficients  for  the  field  reflected  back  into  (B)  are 
the  elements  of  the  first  column  of  the  scattering  matrix  SBB.  The 
field  at  the  aperture  may  be  readily  computed  from  the  modal  expansion 
at  z  =  0,  c  <  x  <  a.  Figure  4.7  illustrates  the  magnitude  of  the 
aperture  field  for  the  array  at  four  angles  of  scan.  The  array 


Figure  4,7  Aperture  field  for  thick-wall  array 
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dimensions  are  a  =  0.62A,  c/a  =  0.12. 

The  method  described  for  the  solution  for  the  thick-wall  array 
is  the  first  combined  application  of  the  MRCT  and  the  generalized 
scattering  matrix  method.  The  principal  advantage  of  the  combined 
technique  is  that  a  geometry  of  some  difficulty  is  reduced  to  two 
geometries  about  which  a  great  deal  is  known.  It  uses  to  advantage 
the  exact  solution  of  the  thin-wall  array  as  well  as  approximate 
solution  to  a  well  known  waveguide  discontinuity  problem.  The 
rapidly  convergent  numerical  technique  for  such  a  combined  geometry 
may  be  quite  useful  for  a  number  of  other  problems  which  can  be 
reduced  to  two  or  more  auxiliary  problems  solvable  by  conventional 
methods  or  by  the  MRCT.  Solutions  for  a  class  of  such  problems, 
specifically  diffracting  surfaces  and  gratings,  are  discussed  in  the 


next  section. 


V.  DIFFRACTING  SURFACES  AND  GRATINGS 


The  modified  residue-calculus  technique  and  the  scattering 
matrix  multiple-reflection  procedure  can  be  combined  to  attack  a 
wide  variety  of  interesting  problems,  such  as  diffraction  by  a  finite 
circular  cylinder,  scattering  from  a  semi- infinite  cylinder  of 
rectangular  cross  section,  and  reflection  at  a  dielectric  step 
discontinuity,  to  cite  only  a  few.  Another  class  of  problems  consists 
of  periodic  diffraction  surfaces  and  gratings.  Since  gratings  are 
finding  current  applications,  for  example  in  resonators  for  lasers 
operating  at  suboptical  frequencies,  alternative  simple  methods  for 
accurately  finding  reflection  and  transmission  properties  of  gratings 
and  surfaces  can  be  quite  useful. 

The  general  form  of  the  periodic  diffracting  surface  treated  in 
this  section  is  shown  in  Figure  5.1.  The  geometry  is  essentially  that 
of  the  array  of  thick  plates  with  a  recessed  dielectric  filling*. 
Because  the  scattering  matrix  representation  of  the  array  of  thick 
plates  is  readily  obtainable  by  the  techniques  of  Section  IV  and 
because  the  scattering  matrices  for  the  junction  between  (B)  and  (B1) 
may  be  found  exactly,  the  combined  geometry  is  most  suitable  for 
description  in  terms  of  the  multiple  reflection  phenomena. 

An  application  of  scattering  matrix  methods  to  diffracting 
surfaces  has  been  made  by  Tseng  [1967],  who  discussed  in  detail 
diffraction  and  surface  wave  propagation  on  infinite  periodic 
corrugated  surfaces.  Although  the  corrugating  fins  were  of  zero 
thickness  in  Tseng's  problem,  the  most  general  unit  cell 
included  multiple  corrugations  of  various  depths. 
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To  avoid  repetition  of  the  work  of  Tseng,  only  the  formal 
analysis  for  a  few  selected  geometries  is  discussed  in  ;his  section, 
primarily  to  indicate  the  extension  of  the  HRCT  and  scattering 
matrix  procedure  to  cases  other  than  the  waveguide  discontinuties  and 
array  problems  treated  previously. 


Analysis 

Let  the  junction  at  z  =  0  in  Figure  5.1  be  designated  by  the 

subscript  1,  so  that  the  scattering  matrices  for  that  junction,  found 

33  3t  ?3  vp 

in  Section  IV,  are  ,  S.,*,  S'  ,  and  S'  .  The  modal  expansion  in 
the  open  region  (F)  is  the  same  as  described  In  Section  IV.  Further, 
let  the  junction  between  (3)  and  (3*)  be  denoted  by  subscript  2.  The 

3B 

scattering  matrix  of  interest  at  that  junction  is  S^  .  If  the 
dielectric  in  (31)  has  relative  dielectric  constant  e^,  the  scattering 
matrix  elements  are  given  by 


(n,m)  =  o®(- 


Z-Z 
n  n« 


B+B 
n  n 


(5.1) 


where 


n 


-(nt/b) 


2 


8n= 


£,  k  2-(m/b)2 

D  O 


are  the  propagation  constants  for  the  mode  in  (3)  and  (B'), 
respectively.  Limiting  cases  of  perfect  electric  and  magnetic 
conductors  in  (B * )  yields 


I 


(5.2) 
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for  electric  wall  at  z=-d,  and 


for  the  magnetic  wall.  The  surface  with  electric  wall  at  z=-d  is 
shown  in  Figure  5.2. 

An  essential  part  of  the  diffracting  surface  is  the  recession 
distance  d.  Propagation  of  the  modes  in  guide  (B)  is  represented  by 
the  propagation  matrix  r,  where 


r(n,m)  =  5^  expC-jB^d) 


(5.4) 


For  example,  propagation  from  z  =  0  to  z=-d,  reflection  at  junction 

BB 

2,  and  propagation  back  from  z=-d  to  z  =  0  is  represented  by  FS2  T. 
Since  the  high-order  modes  are  evanescent,  the  effect  of  the  propaga¬ 
tion  matrix  T  is  to  suitably  weight  them.  The  result  is  that 
truncation  at  small  matrix  sizes  (5x5,  say)  causes  negligible  error. 

The  multiple  reflections  for  a  wave  incident  from  region  (F) 
yield  the  following  expression  for  the  composite  scattering  matrix 


sFF  -  sFF  +  sFBrsf  rsf 


t  sFBrsBBrsf  rsf  rsf 


+  sfrsBBr(sBBrsBBr)2sf  +  ... 


=  sFF  +  sFBrsBBr(i-sBBrsBBr)~1sBF 


(5.5) 
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X 

i 


Figure  5.2  Metallic  modulated  diffracting  surface 


FF 

The  scattering  matrix  S  is  of  course  a  function  of  the  incident 
angle  0,  because  every  scattering  matrix  with  subscript  1  is  an 
implicit  function  of  angle. 

For  the  surface  with  dielectric  filling,  the  transmission  into 
the  dielectric  region  may  also  be  found  by  the  scattering  matrix 
method,  so  that 


sB'r  -  sB’Br(i-sf rsfr^sf 


(5.6) 


An  alternative  interpretation  of  the  geometry  of  Figure  5.1  is  a 
phased  array  with  recessed  dielectric  filling,  so  that  reflection 
coefficients  for  periodic  excitation  from  region  (B1)  may  also  be 
found  by  the  generalized  scattering  matrix  method. 

Strip  grating.  The  diffraction  grating  formed  by  an  infinite 
set  of  thin  parallel  strips  as  shown  in  Figure  5.3  may  also  be 
treated  by  the  multiple-reflection  scattering  matrix  procedure. 

As  with  the  waveguide  diaphragm  of  Section  III,  the  grating  may  be 
reduced  to  two  solvable  auxiliary  geometries  by  considering  symmetric 
and  antisymmetric  excitation  on  either  side  of  the  grating.  Suppose 
the  grating  is  excited  equally  and  symmetrically  from  both  sides. 

For  this  excitation,  the  transverse  magnetic  field  H  in  the  apertures 
of  the  gratings  is  exactly  zero,  so  that  a  magnetic  wall  may  be 
placed  in  each  aperture  without  effect.  As  viewed  from  region  (F), 
the  structure  for  symmetric  excitation  becomes  alternate  magnetic 
and  electric  strips,  which  is  a  limiting  case  of  Figure  5.1  as  the 
dielectric  in  (B1)  becomes  a  magnetic  wall  and  distance  d  reduces  to 
zero.  For  equal  but  antisymmetric  excitation,  the  electric  fields 
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t 


Figure  5.3  Infinite  grating  of  thin  strips 
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cancel  exactly  in  the  aperture  and  the  grating  may  be  replaced  by 

a  perfect  electric  conducting  wall. 

The  solution  for  the  magnetic -electric  strip  grating  may  be 

BB 

obtained  from  equation  (5.5).  However,  the  scattering  matrix  =  I 
(as  in  (5.3))  and  for  d=0,  the  propagation  matrix  also  becomes  T=I. 

In  that  instance,  then. 


(5.7) 


where  the  subscript  s  denotes  symmetric  excitation.  By  inspection, 
the  scattering  matrix  for  the  electric  wall  is 


(5.8) 


where  the  subscript  a  denotes  antisymmetric  excitation. 

For  simultaneous  symmetric  and  antisymmetric  excitation,  the 
incident  waves  from  (E' )  cancel  and  the  waves  from  (F)  add.  As  a 
result,  the  composite  scattering  matrix  for  the  strip  grating  may  be 
written 


KF  *  <F> 


=  |(SFF-I)  (5.9) 

Numerical  results  for  selected  diffracting  surfaces  were  computed 
to  test  the  multiple-reflection  method.  In  all  of  the  test  cases, 
the  reflection  coefficients  were  consistent  with  known  power  require¬ 
ments,  even  with  more  than  one  space  mode  propagating.  Reflection 
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coefficients  for  two  cases  are  given  in  Table  5.1.  The  geometry  is 
that  of  Figure  5.2.  For  both  cases,  the  width  is  a  =  0.75X  and 
depth  d  =  0.5X.  Case  1  has  c=0  and  Case  2  has  c=0.3a.  Case  1  is 
the  geometry  of  Tseng  [19673. 


TABLE  5.1 

Reflection  Coefficient  for  Diffracting  Surface 
(a  =  0.75X,  d  =  0.5X) 


Incident 

Reflection 

Coefficient 

Angle  0 

Case  1. 

c/a  =  0 

Case  2.  c/a  =  0.3 

0° 

0.3845  - 

j0.9231 

0.0177  +  j0.9998 

15° 

0.1586  - 

j0.9873 

0.0732  +  j0.9981 

30° 

-0.2293  - 

j 0.3921 

0.0209  +  jO. 6055 

45° 

-0.3139  - 

j 0.0413 

-0.2406  +  jO. 4956 

Although  the  discussion  of  diffraction  surfaces  in  this  section 
is  rather  brief  and  restricted  to  a  few  basic  geometries,  the  combined 
MRCT  and  scattering  matrix  procedures  may  be  applied  to  a  wide  class 
of  surface  and  grating  problems.  Many  variants  on  the  geometry  of 
Figure  5.1  are  possible,  including  geometries  in  which  reflections 
at  more  than  two  junctions  are  considered. 
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VI.  CONCLUSIONS 

The  analytical  methods  described  in  this  thesis  provide  a  very 
accurate  and  very  rapid  numerical  means  for  solving  a  class  of 
boundary  value  problems  related  to  basic  Wiener-Hopf  geometries. 

As  compared  to  some  variational  or  integral-equation  techniques  for 
solving  the  same  types  of  problems,  the  MRCT  and  the  scattering 
matrix  method,  or  both  together,  yield  the  most  often  required 
information,  i.e.,  scattering  coefficients,  in  a  very  simple  form, 
without  need  for  the  iterative  procedures,  numerical  integrations, 
or  large  order  matrix  inversions  which  often  are  time  consuming 
and  potentially  inaccurate  from  a  computational  standpoint.  A 
particular  advantage  of  the  methods  described  here  is  that  in  a 
modified  Wiener-Hopf  geometry,  a  great  deal  of  information  already 
known  about  the  related  auxiliary  geometry  is  utilized  in  the  more 
complicated  problem.  For  example,  independent  solutions  of  the 
waveguide  step  discontinuity  and  of  the  array  of  thin  plates  may 
readily  be  combined  to  solve  a  much  more  difficult  problem  of  great 
practical  interest,  namely  the  thick-wall  phased  array. 

One  additional  property  of  the  MRCT  is  particularly  significant. 
The  MRCT  automatically  satisfies  the  edge  condition  of  the  problem 
being  solved,  even  when  the  edge  condition  differs  from  that  of  the 
auxiliary  geometry,  because  a  priori  information  about  the  edge 
condition  is  incorporated  into  the  construction  of  the  required 
meromorphic  function.  In  contrast,  satisfaction  of  the  edge 
condition  in  solution  of  a  doubly  infinite  set  of  equations  by 


truncation  is  strongly  dependent  upon  the  truncation  procedure 
(Mittra  [1963]). 

Although  a  variety  of  waveguide  discontinuities  were  solved  by 
the  MRCT  ,  combination  with  the  scattering  matrix  procedure  permits 
solution  of  such  problems  as  thick  diaphragms,  trifurcations  (Pace 
[1964]),  multiple  diaphragms,  and  rectangular  posts  in  waveguides. 
Also,  the  MRCT  can  be  extended  to  circular  cylindrical  geometries, 
so  that  a  variety  of  corresponding  discontinuities  in  circular 
waveguides  .may  be  solved. 

The  applicability  of  the  analytical  methods  of  this  thesis  to 
a  wide  variety  of  problems  in  both  open  and  closed  regions  will 
permit  numerical  solutions  in  a  convenient  form  for  many  problems  of 


current  interest. 
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